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In case of configurations of a cluster that are mixed, 

the question of interactions between two different kinds of 

bonds arises. For example in case of the Ga^InT As alloys, 

As ^ 

a configuration of the kind may be a part of a 

Qsy ^'In 

cluster. The interaction between these bonds 1 and 2 

A 

is of a kind but between a Ga-As and an In- As bond. It is 
demonstrated in the following that lliis interaction is the 
average of H^‘s corresponding to pure Ga-As end pure In-As 
crystals, as a reasonable approximation. 

Let bonds 1- and 2 be described by quantum states }b^> 
and jb 2 > respectively. Then, 

|b^> = (jAs^> + lGa^> ) -|- 

and |b 2 > = ( |As 2 > + |ln 2 > ) 

in terms of the sp^ hybridised orbitals at the various sites 
directed along the respective bonds and N is the nomialisa- 
tion constant ~ 2, assimaed to be the same for the Ga-As and 

In-As bonds. We are concerned with , i.e. 

<1^1 H {iAs2> + |ln2>}] 

= f^As-j^ jH |A s 2^ + <Ga2^lHjAs2^ 

+ <As3^|H|ln2> + <Ga^jH|ln2>3 
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probabilities are indicated in Fig. 3*7(c). Configurations 
and are identical but are listed separately to 

retain symmetry imder exchange of the alloyed ions Ga and In, 
Now, 



where and G^ are the matrix elements of the CPA Green’ s j 

operator. are similarly determined and the n x n matrix ; 

j^exact for each configuration by substituting in the corre4 

parameters depending on vdiether the bonds i and j are the Ga— As j 

or the In-As type. An average over the trace of the averaged 
Green's operator matrix is taken to yield the averaged diagonal 
Green’s operator sought, since the four bonds are not symmetrical 
in all the configurations. 

(d) Next, the seven bond cluster consisting of a bond and 
its six nearest neighbours is considered. The matrices involved 
are therefore 7 x' 7. Eight configurations and their respective 
statistical probabilities of occurrence are shown in Fig. 3.7(d). 
Matrix elements of the CPA Green's operator between common cation 
and common anion adjacent bonds, second nearest neighbouring 
parallel and non-parallel bonds, and the diagonal elements all 
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CHAPTER IV 


THE CHEMICAL PSEUDOPOTENTIAL THEORY, APPLICATION 

AND CCMPARISONS 


4.1 Introduction ; 

The parameters used to describe the III-V semicon- 
ductors in Chapter III have been obtained by fitting with 
experimental observations of the photoelectric thresholds 
and x-ray photomission spectra. Tlie procedure is unsatis- 
factory due to two main reasons. Firstly, experimental 
results obtained by different workers differ. So much so that 
differences in parameters for the same alloy obtained from 
different sources vary more than parameters corresponding to 
different alloys (See Table Chen and Sher (1978) in fact, 

first determine parameters corresponding to the six III-V 
semiconducting materials considered, separately from two or 
three different sources of experimental observations differing 
sometimes widely. Finally, the best fitted values are 
determined by a hit-and-trial method, such that after an 
energy dependent broadening is imposed, the density of states 
function obtained best resembles the esqperimental x-ray photo- 
emission spectra curves. This involves much airbitrariness. 
Secondly, the bond basis set is non-orthogonal, with conside- 
rable overlap between bonds on different pairs of atoms. 
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where T is the kinetic energy operator and V is tlie effective 

€LC 

bond potential. A functional form for the V may he obtained 
in terms of the atomic potentials and V at the sites 
labelled by a and c, and the corresponding orbital functions 
and 0^ respectively which constitute the bond under con- 
sideration, by writing Eq. (4,14) as 

(T + [(a> + |c>] « C|a> + jc>] 

the normalisation constant being irrelevant here. This can be 
simplified to 


T^|a> + T^|c> + Vg^^[|a> + |c>] = e^[ja> + (c>] 


i.e. na> + |c>3 = 6^[|a>+{c>3 


c c 


ac 


\c " + 




(4.15) 


where T_ and are the kinetic energy operators corresponding 
to the electrons bound to the atoms at sites a and c respec- 
tively. In case of carbon and silicon in a diamond lattice 
this potential along the bond is symmetric (Bullett 1975) 
while in case of III-V semiconductors forming a zinc-blende 
structure, the potential is asymmetric. 


The pseudopotential off diagonal matrix elements are 

therefore like <a’ c’ fV^®^^}ac>, where ja’c’> is a bond orbital 

a’ c’ 

nei^bouring to |ac>. To evaluate these, one could use the 
’ Wigner trick ' similar to the one used by Anderson for his 
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SYNOPSIS 



The study of metallic alloys like brass, steels and 
white metal has had a long history. Semi conductors and 
their alloys like Si, Ge, GaAs, Iniis and GaAs^ ^1-x 
attracted attention with the advent of photovoltaic cells, 
as have amorphous semi-conductors like hydrogenated a-Si 
films as proposed cheap, but efficient solar cells. For a 
long time the study has been empirical, essentially a metal- 
lurgical approach. Till recently there has been a lack of 
xmderlying imity of theory, and consequently controversy 
over seemingly contradictory experimental data. 

Often, the properties of the alloys and the consti- 
tuents or that of the crystalline and the amorphous materials 
differ drastically. In order to imder stand why these systems 
exhibit these properties and to be able to predict specific 
properties, not only must one understand how these arise in 
the pure or ordered systems but also how they change on 
introducing disorder. Clearly, this requires a microscopic 
theoretical study that is systematic. The thesis is an 
attempt in this direction. CHir attention is confined to 
the electronic properties of disordered systems. 

A self consistent cluster CPA formalism to obtain a 
configurationally averaged Green' s function for electronic 



( 11 ) 


Fig. 4,9t (1) -In 0^*^^ and (11) -Im 0®*** for 

GaQ^^ ^^0.5 ^*1^ ehaalcal bond 

paeudopotantlal paraaetara. 
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tree is always more diffusive than the same on the real 
lattice from which the tree was formed by ' delinking' • 
Estimates of the localisation length in the CTA are thus 
always over-estimates, as are the estimates of the mobility 
edges* 

We can rewrite Eqn, (5*1) Bomevfliat. Following Econcwnou 
and Cohen (1972), 

OqCz) being the self energy, 

then, 

if* lOkJriC® In) 

» / dE lim {nQ(E)/[l + — ]} 

"" n ”*■0’*' 

We have used the herglotz property of the Green’s 
function and the self energy. 


i*e* sgn Gqq(Z) 3 « •• sgn [Im Z)| 
lim Im Qqq (E - in) = n^CE) 

n **0+ 


the local density of states; and since Oq(Z) is analytic 
in the entire upper half plane In Z > 0, 

Im cf Q(E-in ) 


lim 

n"*'0^ 


da q(E) 

- — jjg 


Now, wiiSiin the band, nQ(E) / 0* If the state is extended, 
d<y^d& diverges and Pq(E) « 0. On the other hand, in the 
localised domain, dcf^dB is finite and *♦ 0 when the degree 
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spectra is developed in C!hapter III, This is based on the 
Augmented Space method of Mookerjee, The method is superior 
to the CPA (reviewed in Chapter II) in that statistical clus- 
tering effects and off-diagonal disorder over finite cltisters 
are successfully taken into account without destroying the 
herglotz property of the Green 's function. This is essential 
for physical relevance. 

The theory is applied to the valence bands of subs- 
titutional alloys, in particular to those consisting of 
pairs of GaAs, InAs, GaSb, InSb, GaP, InP with a common 

0 ^ X £ 1. These materials form tetrahedrally coordinated, 
strong covalently bonded solids in a zinc-blende lattice 
structure. In the alloys the fee sublattice of the common 
ion is ’ regular • , while that of the other ion is randomly 
occupied by the constituent atoms. 

These are studied in the Bond Orbital Model ( BOl^ 
which has the advantage over the Site Orbital Model of having 
rapidly falling off-diagonal terms. In doing this we are 
neglecting the bonding-antibonding overlaps. Also the anti- 
bonding orbital bases have been found to be inappropriate 
for the study of the conduction band, and we do not get any 
estimates of the band gaps. A parameterised hamiltonian 
fitted to experimental data is used to obtain the Green's 


cation or anion mixed in various properties e.g. Ggi I 




+ ?:nuunR«ll»eJ Prob«hiUty 
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functions for the ternary alloys in the Virtual Crystal 
Approximation (VGA). A Ray Diagram technique, described 
in Appendix A, is employed to evaluate density of states 
and consequently the complete Green's functions by the use 
of Kramers-Kronig relations. A similar technique is used to ; 
obtain the off-diagonal elements of the Green’s operator, 
required for the cluster embedding in the effective medium. 

The VGA result is used as a starting point for the subse- 
quent CPA calculation. Cluster CPA Green's functions corres-; 

I 

ponding to exact clusters of various kinds embedded in the | 

I 

CP-medium are then obtained. A weighted average of these j 
gives the required configurationally averaged Green’s functioi 

I 

Results are obtained for one-bond, four-bond, seven-bond and ■ 

I 

six-bond-ring clusters. The CP-medium is not made selP-consii 
tent with respect to the clusters larger than one-bond. 1 

Consequently, there is no improvement of the band edges or th< 
structure near them. But the cluster characteristics near thj 
band centres are reproduced. In addition, /hypothetical alloyij 

. . j 

with widely differing potentials are studied, revealing the j 
'impurity bands ' . These show considerable sensitivity to 1 

I 

clustering and off-diagonal disorder. 

The parameter fitting procedure to fix the hamiltonian 
carried out using the photoelectric thresholds and x-ray 
photoemission spectra involves a degree of arbitrariness 
depending on different experimental data, often differing 
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shown to effectively traiisform the Hill determinant such 
that, for every eigen energy, most of the contribution is 
lumped in one term on the diagonal. Opening the Hill deter- 
minant about this term gives exceedingly fast converging 
perturbation expansion for the eigenvalues and eigenfunctions 
The boundary layer in the (k,?^ ) plane which has been 
believed to be unsurpassable by all other methods available 
before, is crossed with ease, for all anharmonicities. Furthei 
more the convergence is so fast and the accuracy limited only ; 
by the machine precision, that the solutions can be called 

’ exact in a numerical sense. ' 

I 

Other problems that can be attempted are a) the double-' 
hump potentials eg. in the ammonia molecule, b) perturbations 1 
as by electric fields i.e. Stack effect, c) anharmonicities othe; 
than the polynomial ones, etc. The approach is very promising 
and will probably prove very useful in the study of other non- : 

i 

linear problems. 
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with each other considerably. The non-orthogonality of the 
bond bases is claimed to be absorbed in the parameters. This 
is an approximation of unclear validity. In CSiapter IV, we 
have developed a chemical pseudopotential scheme to determine 
the parameters, starting from Herman-Skillman tables of 
self-consistent Hartree-Fock one-electron potentials and 
orbitals. The chemical pseudopotential method not only takes 
the non-orthogonality of the bases into accomt but is illus- 
trative of the chemistry of the bonding in a more transparent 
fashion. Furthermore, the ab initio procedure is independent 
of parameter fitting to experiments. The VGA, CPA and the 
cluster-CPA calculations are likewise made on ttie same systems 
and comparisons made. 

In both cases, the cluster CPA alters essentially the 
top of the valence band. This is to be expected because the 
inner states must see less of the details of the neighbourhood. 

In Chapter V, the question of the nature of electronic 
states in disordered systems is tackled. A criterion for 
localisation involving the convergence or divergence of the 
averaged localisation length is introduced. An estimate of 
localisation lengths of electron states in weakly disordered 
Cayley trees is made. Results are exact near the mobility 
edge, producing the linear dependence on energy expected in 
effective medium theories. In general, nonlinear integral 



Eqn. (A.1) for each TT in each tetrahedron of the IBZ and 
accijmulate the sum, (vii) Use Kramers-Kronig relations to evali- 

I 

the complete function F(E). 

VJhen f^) in the integrand does not reflect the symme- 
tries of the BZ, contributions in each of the other parts of 
BZ similar to the IBZ, obtainable by the application of a symme 
transformation of the complete BZ, may be obtained by inserting 
the value of f (k) at points corresponding to the in the IBZ 
scanned. Since. the Hamiltonian is necessarily symmetric imder 
these very transformations, all other factors remain the same.. 

A. 3 G eneralisation and Application to Off-I)Lagonal Green's 
Functions : 

The BZ of an fee lattice is shown in Flg.A.3in which 
so6ie symmetry points and directions are also indicated. Label^ 
in the band structures for GaAs and InAs on Fig. (3.5) corres- 
pond to those shown here. The IBZ is shown in FigJl,4'^lch 
can conveniently be divided into three tetrahedrfl 1,2 and 3» 
all apexed at r » the centre of the BZ (Fig. A. 5)* coordi- 
nates shown are in the units of ii/l6a, where a is the edge of 
the direct fee lattice cubic omit cell, t corresponds to k = 0 

and i s.theref or e , taken to be the origin of "ttie coordinate 

^ \ 

I 

system, and the coordinate axes X,y,z pass through the centres} 
of the square faces of the BZ (Fig. A.4)* Point A and the 
symmetry points r , L, U, X, K are the vertices of the IBZ, 

LUAK is a sixth of the hexagonal face, while UXA is an eighth 
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equations are obtained. Close to the band edge, solutions 
necessitated numerical computation. This has been carried 
out and the results presented. 

As a case of divergent but asymptotic perturbation, 
the problem of the anharmonic oscillator is also considered. 

An elegant non-perturbative method using the 'difference equa- 
tion, to evaluate eigenvalues and eigenfunctions is presented. 
The work, done in collaboration with Er. K. Banerjee, S.S.Kanwal 
and S,P, Bhatnagar [Proc. Roy. Soc,, London, A g60. 575-586 

(1978)] does not get worse with higher quantiom numbers for, 

2 A- 

sa^r, the quartic anharmonicity, i.e., V = kx +Xx , or with 
the anharmonicit 3 J’ parameter X*The hitherto imsurpassable 
’boundary layer’ in the (k, x) plane is crossed with ease, The 
calculation takes considerably less time than the best methods 
otherwise available only for some regimes of (k, x) and few 
quantum numbers. Fast quadratic convergence and accxiracy 
limited only by the ’PRECISION* of the conqjuter calculation 
are other highlights of the formalian. Some eigenvalues accu- 
rate to 15 significant figures are tabulated, calculated in 
the Double Precision mode on IBM 7044 and EEC 10, 



CHAPTER I 


INTRODUCTION 


1,1 Disordered Systems s 

Until recently, most of the attention of physicists 
has been occupied by regular or ordered systems. Yet there 
is no dearth of those essentially disordered. Examples 
abound. Besides the geophysical and astrophysical systems, 
liquids, amorphous solids, glasses, random alloys, polymers 
are seme of them. 

The original metallurgical approach for discovering 
useful metal alloys, for example, was to ' mix and see ’ if a 
useful, reproducible property is obtained. In this way 
steels, the discovery of which revolutionised Industry, were 
obtained by alloying iron with small amounts of carbon, 
chromium, etc. Varied proportions of the same constituents 
gave steels exhibiting various degrees of hardness, and 
other associated properties like ductilit 3 /- and tenacity. 
Similarly, brasses, alloys of copper and zinc, have been 
in use. White metals too are verjr useful as antifriction 
coatings for bearings. Of these, the Babbit metals are composed 
of mainly tin alloyed with smaller amounts of antimony and 
copper, while the Brittania metals also contain zinc, lead 
and bismuth. The latter are also used for maJfclng table 


ware. 
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K. Banerjee and others 

eigenfunction Just begins to increase in magnitude. The contribulacMni to the moments 
firom the rest of the configuration space is estimated to be ies than 0 { 10-“) fmr all 

the mmnenta of the table 3 ; (in) the moments may now bo ezjneamd as sums of 
n 

integrals of the type a:**e-“**da:, {k ■» integer), which are otouned reouravely 

starting from the incomplete Gaussian integral J*e-***dar. The computation of 

moments therefore requires no integration; {tv) each moment in table 3 has been 
checked by varying a in the range (7), 

The matnx elements of the operators z and ®* in the lowest ten states of the 
quartic oscillator {H » |»»+®<) are given in table 3. The best known ^imates 
(Chan As Stelman 1963; Reid 1970) of the transition moments of the quartic 
osciUator were obtained by using variational eigenfunctions. The variational 
eigenfunctions ari known to be much less accurate than the corresponding eigen- 
values and are unsuitable for the computation of moments. Thus only the first 
few (two or three) non-aero moments for the transitions wme obtained in these 
works to some reliability (6 figures). We emphasize that ail moments except the 
lowest one for any transition may be obtained from (12), without integration. 


6. Dtsonssioir 

(i) The confidence in the accuracy of the computed eigenvalues is derived from 
the following checks: (a) The computed eigenvalues are stable with respect to 
changes in the scaling a (see ( 7 )) and for different initial estimates for the energies; 
(6) Three separate computations with' an increasit^, a decreasing, and a nearly 
flat, sequence of determinants (see the discussion after (11)) yielded eigenvalure 
differing by at most two units in the 16 th significant figure; (e) The upper and lower 
bounds for the eigenvalues differing in the 16 th figmre by at most two units were 
obtained by the ‘sign check’ described in §2. 

(ii) It may be emphasized that no integration or diagonaiizalion is imcessary in 
this method which makes it attractive for the eigenvalue problems of tbe Uneiff 
operators. 

(iii) In perturbation theory with A as the small parameter, this problem beltmgs 
to the singular perturbation class. The uniform applicability of this method for all 
A underUnes its non-perturbativo character. For this reason the method is expected 
to apply to i«t»blems involving large perturbations. 

(iv) The need for introducing a scaled basis was resized by Reid (1970) udio 
used the linear variation method with harmonic oscillator basis for the quartic 
oscillator eigenvalue problem. However, in a varfarional framework tl» use of a 
scaled basis becomes introctably laborious for the following reason. In a variational 
cominitation the first n (say) eigenvalues are obtained together. Since the appro- 
priate scaling is different for different n and A a single scaling is not suitable fior the 
computation of all n eigenvalues. A compromise scaling must therefore be used. 





Semiconducting alloys e.g. GaAs-GaAlAs 

have recently become Important because of their use in solar 
cells in search of an alternative energy source. Also 
GaAsP cells are used as light emitters and HgCdTe as 
IF-detectors. 

The special properties of each of the alloys mentioned 
above often differ drastically from those of the constituents. 
In order to really understand why the alloys esdiibit these 
properties and to be able to predict specific properties, 
not only must one know how these properties arise in the 
pure constituents, but also how •ttiese change on alloying. 

For example, in case of solar cell materials, how the 
band gap changes on alloying is an important factor. Of 
co^lrse, in general, the shapes of the bands together with the 
band gap, determine the optical and optico-electrical proper- 
ties. The knowledge of how these change on alloying would 
make it possible to determine the optimimi proportions of 
constituents in an alloy, that would exhibit the desired 
properties and behaviour, in a theoretically systematic 
manner, instead of the ' hit and trial' experimental method. 
Clearly, this reqxiires a study at the microscopic level. 

This is the view point of the theoretical physicists. 

Interest in amorphous states of pure and mixed intrinsic 
semi-conductors eg a-Si and Te^ ^®10 sparked 

off by Ovshinsky's (1968) reported discovery of its switching 



property, A-Si now shows tremendous promise in construc- 
tion of cheap and efficient solar cells. Although solar 
cells using crystalline Si have been designed fairly success- 
fully in the laboratory, these are not commercially viable, as 
tedious methods of preparation of large enough crystalline cells 
are prohibitively expensive. Solar cells using a-Si, on the 
other hand, are comparatively non-expensive both because of 
simpler methods of preparation and because thin films require 
very little material. Thus large areas coated with these 
films may be easily obtained. Such cells v^ith a wor king 
efficiency of 4^ are now available commercially. Those 
working at about 8/3 efficiency have also been reported but 
have short lives. Besides, these are not easily reproduced. 

To understand these features, again a microscopic analysis 
is indicated. 

Glasses, too, are disordered materials in the amorphous 
state. Other systems like amorphous Sulphur, Selenium, 
Tellurium,etc, have also been extensively experimented with, 
though a microscopic understanding of the behaviour of these 
is yet incomplete. 

E5speriments to detect presence of voids and dang i i ng 
bonds, sensitivity to annealing, judged by measuring changes 
in the ESR response, conductivity, voids, stresses, diamag- 
netism, etc, as the annealing rates are varied, seem to indi- 
cate inherent differences between some materials. Some 



materials defy preparation by the method of quenching the 
liquid state, always ending up in polycrystalline states 
this way, while some others should better be called ' super 
cooled liquids’ because no latent heat aroxmd the apparent 
change of state is involved, Puirther, some materials do, 
while others do not even show a ’ glass transition’ • Some- 
times the essential difference between some materials occurs 
because some materials can be made in bulk, while some are 
obtained only as films, ihe behavioiir of the former is 
governed by bulk properties, while that of latter must 
essentially be due to surface properties. In fact, studies 
on silicon-dioxide films have revealed that some properties 
of silicon samples are truly due to the silicon-dioxide layer 
that gets very easil3'‘ formed on the sirrface due to oxidisation. 
Surface properties are, likewise important for transition 
metal alloys in their use as catalysts for chemical reactions. 
Spin glasses, met-glasses, garnets and other random magnetic 
systems are some of the other systems of concern. 

Such varied behaviour makes a xmiversal theory for all 
disordered systems virtually impossible. Those mentioned 
above, and others studied experimentally with a view to reach 
a better understanding at large, therefore, need to be 
classified. 



1.2 Classification? 


1) Compositionallv Eisordered systems, i.e. those in which 
the locations of sites fonn a lattice, occupied randomly by 
atoms of the constituents. These include (a) defects and 
impxirities in otherwise periodic sj^-stems, e.g. coloiar centres 
in ionic alloys, (b) metal alloys i.e. the brasses, steels, 
white metals, and (c) the semi-conducting alloys, e.g. 

2) structurally disordered systems include the amorphous 
materials. These ares 

a) Dielectric films like Si ©2 or ionic insulators like 
^ 2 ° 5 * 

b) The amorphous semiconductors, which are further sub- 

divided into two classes - (i) the tetrahedrallj? 
bonded (TB) type e.g. a-Si, a-Ge, a-GaAs, and (ii) the 
lone pair (LP) type, e.g. a-Se, a.-~ks2^2* sire 

distinguished because of general differences shown 
experimentally. For example, the former are generally 
found only as thin films containing voids and dangling 
bonds, while the latter, formed in bulk, do not have 
enough density of these to be detected bjr the porosity 
and ESR experiments respectively. Again, the TB 
a-semiconductors are very sensitive to annealing 
indicated by the disappearance of the ESR response, 
increase in resistivity, decrease in the number of voids 
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and relaxation of internal stress. These also show 
considerable diamagnetic enhancement in the amorphous 
state. !The LP type of a-semiconductors, on the other 
hand, are much less sensitive to annealing, show no 
stress thereafter, and the diamagnetic enhancement 
is not much. Finally, these exhibit. a glass transition, 
while the TB type can be obtained only as films and 
show no glass transition, passing into the crystalline 
state instead, 

c) Quasi one-dimensional systems like TTN-TCNQ have been 
designed and dealt with, essentially to gain under- 
standing, since theoretically the problem is relatively?’ 
simple. 

3) Polymers are a subclass of structural disorder because 
of their structural peculiarities. These consist of long 
linear molecules which fill space. Interactions between 
various parts of a chain as well as other chains exist. These 
often exhibit a gelling phenomena related to the localisation 
in disordered systems, vdiich must be xmderstood. Interest 

in polymers arises from analogous systems like plastics in 
material science and’ gels, polymer solutions, Dl^, the 
genetic material etc, in life sciences. 

4) Magnetic disorder occurs when the spin orientation is 
random forming random para- and ferro-magnets. Five systems. 



viz. Cu-Mn, Au-Fe, Au-Cr, Ag-Cr and Ag-Mn are crystalline 
random alloys consisting of moderate concentrations of 
the magnetic constituent and show a new magnetic phase 
called the spin glass • Kondo effect, i.e. disappearance 
of magnetisation below a certain temperature, of very 
dilute magnetic impurities in a non-magnetic host, is an 
unusual property esdiibited by this class. Finally, met- 
glasses, garnets etc. consist of a structurally disordered 
matrix of transition metal atoms embedded in a glass-forming 
matrix, e.g. 

1.3 No n-Reproducibility of Experimental Result ^s 

Observations made on the same system by different 
experimentalists sometimes do not tally. Of course, a 
particular sample can accidentally develop a specific 
anomaly during its preparation or while observations are 
made on it, therefore, an experimentalist takes the pre- 
caution of repeating observations on a sample, as well as 
various samples prepared under the same conditions before 
quoting the results. 

Some parameters that affect the properties of the 
sample obtained are (l) method of preparation, (2) the 
environment, (3) the source of raw material and (4) the 
consequent analysis of data. 

For example, while lone pair amorphous semiconductors 
can be prepared in bulk by quenching the liquid state. 
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tetrahedrally bonded amorphous semiconductors fona poly- 
crystalline samples in this process, exhibiting very different 
properties. Again, thin film of a-Si obtained from rf- 
sputtering in an inert gas environment is now distinguished 
from that obtained by glow discharge through silane (SiH^) 
gas. The former contains microvoids, with dangling bonds 
in them, whereas the latter, referred to as h 3 rdrogenated 
a-Si, is different essentially because the available H-atoms 
enter the sample and satisfy the dangling bonds. Presence 
of impiirities in dilute amounts also can alter properties 
of a sample, hence the source of raw material used and the 
environment of preparation must be taken into accoimt while 
trying to imderstand the mechanism of behaviour. Finally, 
the analysis employed to interpret the observations can intro- 
duce controversies. For ex:ample, band gaps as determined 
from electrical and optical properties need not give similar 
results. 

There is, therefore, an •urgent need to correctly inter- 
pret experimental results and clearly understand the various 
parameters which affect their outccane, in a systematically 
developed theoretical approach. 

1.4 Theoretical View ? 

The solution of the many body problem of structures 
and properties of crystalline materials is facilitated to 
a great extent by the Bloch’s theorem. Lattice periodicity. 
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or translational invariance of the potential is central to 
the theorem. As a consequence, in the momentum space, the 
degrees of freedom reduce to a few, and the problem becomes 
tractable. 

The idea of bands and band gaps came out of this. Much 
of solid state phenomena of crystals is understood in terms 
of these. For example, the position of the Fermi level 
within the bands or gaps determines the conductive nature 
of the material. For the lying within the band gap AE, 
if aE ~ kT, the material is semiconducting. For aE « kT or 
overlapping bands, metallic behaviour is produced. For 
aE » kT, we have insulators, which are transparent if aE 
is greater than the energy of visible photons. But glass, an 
amorphous material lacking lattice s^nnmetrj'-, must have a 
large band gap, because of its transparency to light. Thus, 
the general belief that Bloch *s theorem or lattice symmetry 
is necessary for the existence of' band gaps, is untrue. 

Actually, the long range order that Bloch’s theorem 
ass'umes, is not essential for the existence of bands and gaps 
between them. That well defined bonding and antibonding 
energy levels in isolated molecules spread into bands in the 
solid state has been the view of chemists. In fact, all 
band structure calctfLations of the density of states take 
accoimt only of interactions between near neighbours. To 
sum up, long range order is not necessary, while the local 
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Fig. 1.1; A typical radial distribution function f(r) 
scaled by r2 for crystalline (full lines) 
and amorphous materials (broken lines), where 
f (r) ; the number of atoms around distance r 
per unit volume.- 

r ? the n-th near neighbour distance in 
crystals, and 

r : the distance within which all n-th 
near neighbour are situated in both 
cases. 
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bonding behaviour is sufficient to determine the bands j 

and properties arising from these. i 

Experiments on amorphous materials reveal that 

I 

the coordination number is retained. Typical radial 
distribution functions f(r), the number of atoms round 

p 

the distance r per unit volme, scaled by r to remove 

I 

the inevitable steady rise in f(r), for crystalline and ; 

amorphous states of materials like Ge, Si are shown in 

Fig, 1,1, In the amorphous states, the peak at r£ , where | 

i 

r^ are the n-th near~neighbour distances, though gets . 

broadened, does not shift. That the curve gets smooth | 

earlier -Uian for crystals, is indicative of the absence : 

j 

of long range order in amorphous materials. This obser- i 

vation supposrts the above view and suggests a local theory 
based on bonding that must be attempted. ! 

In general, one would like to look into the nature i 

I 

of various excitations, e,g. electrons, phonons, magnons, { 

I 

etc. As Elliott et.al. (1974) have shown, within the ; 

I 

Green* s operator formalism, studies of each of these can be i 

j 

cast on very similar frameworks. Methods developed for one are, 
thereafter, easily extended for other excitations . In the followini 
work we concentrate on the electronic properties, in particular, 

i 
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the density of states. In Chapter II, we review earlier 
attempts at developing the theory of electronic structure 
of random alloys. The idea is to examine in detail the 
drawbacks of these approaches so as to form a basis for our 
generalisation. The Coherent Potential Approximation (CPA), 
which has been the most successful developnent so far, has 
been examined with a view of determining the best avenue open 
for generalisation. Unsuccessful or pairfcially successful 
attempts at generalisation have also been discussed. Finally, 
an analytically sound cluster CP - generalisation has been 
developed, which at the same time is tractable for reason- 
ably realistic calculations on real alloys. This generali- 
sation is based on the augmented space approach introduced 
by Mookerjee (1973) and later developed by various authors 
(Kaplan, Gray 1976, 1977, Kaplan, Leath,Gray, Diehl I960), 

This is presented in Chapter III and applied to III-V semi- 
conducting ternary alloys obtained from pairs of GaAs, InAs, 
GaSb, InSb, GaP and InP, with a common cation or an anion in 
various ratios. These are covalently bonded with a four-fold 
coordination, and fall in the ‘category of compositionally dis- 
ordered systems, i.e, have a randomly occupied lattice. The 
approach is basically a tight-binding formalism suitable for 

our theory. Bond orbitals built out of symmetric combinations 
3 

of sp^-hybridised orbitals at each atom site are used as the 
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bases. Matrix elements of the Hamiltonian are used as 
parameters. These are determined from fitting to specific 
experimental results. 

In Chapter IV, a chemical pseudopotential method for 
an ab initio derivation of the crystal Hamiltonian is deve- 
loped. This is a generalisation of Bullett's (1975) appli- 
cation of the pseudopotential theory to bond orbitals of 
systems like crystalline carbon and silicon. The method 
generates the spectrum from atomic potentials, takes into 
account the non-orthogonality of the basis and leaves us 
free to predict experimental results, rather than fit to 
them. The spectra corresponding to the pure constituents 
are then used as input in the cluster-CPA calculation for the 
alloy spectra. 

The question of the nature of the electron states in 
disordered systems, diffusive or localised, pioneered by 
Anderson (1958) has been tackled in Chapter V. Anderson 
suggested that localised or non-diffusive electron states may 
exist in disordered systems. Much confusion prevailed regard- 
ing the existence of such states. However, presently, there 
is general acceptance of the existence of localised states 
and a bulk of experimental data may be explained on the basis 
of this assumption. Borland (1964) and Ishii (1973) showed 
that in one-dimension, all states become localised for any 
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finite disorder. Mott (1967,68.) introduced the idea of sharp 
inobilit3?' edges separating extended and localised states in 
three dimensions and used this concept to explain the elec- 
trical properties of amorphous semiconductors. As the degree 
of disorder increases, the mobility’- edges move into the band. 

At a critical disorder all states become localised. This is 
the *inderson transition. 

We develop a criterion for the calculation of the 
averaged localisation length of electron states and estimate 
these in weakly disordered systems in the Cayley Tree Approxi- 
mation, The problem is exactly solvable near the mobility 
edge giving the linear dependence of the averaged localisation 
length on the energy as e 3 q)ected in effective medium theories. 
Away from the mobility edges, the localisation length has 
been determined numerically. 

As a case of divergent but as3niiptotic perturbation, 
the problem of the anharmonic oscillator is also considered, 
iin elegant non-perturbative method using the difference equa- 
tion, to evaluate eigenvalues and eigenfmctions is presented. 
The work, done in collaboration with Dr. K, Banerjee, S.S.Kanwal 

and S.P. Bhatnagar, does not get worse with higher quantum 

P A 

nijmbers for, say, the quartic anharmonicity, i.e., V = kx + Xx , 
or with the anharmonicity parameter X, The hitherto unsiirpass- 
able boundary layer in the (K, X ) plane is crossed with ease. 
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The calculation takes considerably less time than the best 
methods otherwise available only for some regimes of (k, \ ) 
and few quantxam numbers. Fast quadratic convergence and 
accuracy limited only b}^ the ‘ PRECISIOi^’ of the computer 
calculation are other highlights bf the formalism. Some 
eigenvalues accurate to 15 significant figures are tabulated, 
calculated in the Double Precision mode on IH'I 7044 and 
EEC 10. 



CHAPTER II 


CRITICAL REVIEW OF E.:iRLiER THEORIES OF ALLOYS 

2.1 Introduction 8 

In this chapter we shall review critically earlier 
attempts to calculate density of electronic states in dis- 
ordered systems. Though attractive at first, some approaches 
become intractable for a realistic calculation, while others 
involving simplifying assumptions either become far removed 
from reality or lose features essential for physical relevence. 
The objective, therefore, is to highlight shortcomings of 
these attempts, while at the same time discuss the basic ideas 
behind the different approaches, so as to expose avenues of 
generalisation, which not only conform to physical reality 
but are tractable for actual application. 

2.1.1 Rigid Band Model (REM) : 

In this model, the density of states function is 
assumed to be uniform in the alloy. The electronic concen- 
tration in the conduction band is altered and set by equalising 
the Fermi energy at each atomic site. Charge transfer therefore, 
must occur, although experiments on Cu-Ni alloys seem to indicate 
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its absence (Wenger et al 1971, Love et al 1971, Bennett et al 
197^). The model predicts disappearance of ferromagnetism as 
concentration of Ni is increased (Kirkpatrick et al 1969) » 
since the partially filled d-bands in pure Ni responsible for 
ferromagnetism would have to be filled when alloyed with Cu 
having one more electron and hence a higher Fermi level. 
Qualitative agreement with experiments is found, . though not 
quantitatively (Mott and Jones 1979)* The RB'I predicts this 
disappearance to occur at Cu concentration of 10^ while it 
actually occurs at 6055, The model, however, agrees well even 
quantitatively with observed optical gaps and manjr features of 
the Slater Pauling curves of satiiration magnetisation vs, 
electron concentration for some of the alloys of transition 
metal series (Hill and Matthias 1968) , 

Needless to saj'', the RBM is too crude for cases when 
distinct bands corresponding to each of the constituents of 
the alloy are present, as when many properties of the alloy 
sensitive to the constituents are of concern. 

2.1.2 V irtual Crystal Approximation (VCA ) : 

In this model, the actual random alloy potential is 
replaced by an average, 

7 = O * (1 - O) Vg 

if and Vg are the potentials of the constituents A and B 
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present in atomic concentrations c and (1-c) respective!}^. 
The model Hamiltonian, 


H = E e . + S V. . 
i ij 

is then translationally invariant and possesses the crystal 
symmetries. Band structure and other properties within this 
model can then be evaluated using standard techniques for 
crystals. 


The VGA cannot be expected to reveal any features 
characteristic of disorder other than the crude averages. 
However, it does serve well as a zeroeth approximation for more 
refined theories. In case the potentials of the constituents 
are very similar, it is not a bad approximation e.g. Ge-Si 
alloys, but in cases like Cu-Zn and Li-Hg alloys, or when 
we are dealing with impurity bands, tally with experimental 
observations is not even satisfactory. The model can be said 
to be valid when the mean free path is large. Valence bands 
of III - V semiconducting alloys studies in this work have been 
dealt within this model in the past. Some features like Van 


Hove singularities persist, which 
alloys (Fig. 3.9). 

Also in general, alloys 
exhibit bowing of the band gap 
on alloying (Fig. 2.1). Neither 
the RBM nor the VGA show this 
bowing behaviour. 


are not expected in the real 



Fig. 2,1: 'Bowing' behaviour 
of the band gap 
alloying. 
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2«1«3 Minimum Polarity Model (MFM) : 

Van Vleck (1953) imposed electrical neutrality at 
each atomic site as an ansatz- for understanding ferromag- 
netism in Ni« Negligible charge transfer observed for Ni-Cu 
alloys indicates applicability of the model. But clearly, the 
model has limited scope while agreement with experimental 
observations is crude, 

2.1.4 Virtual Bound State Model (VBSM) : 

Also called the Friedel-Anderson model, this really is 
between the extremes of the REM and the MEM and is applicable 
for transition metal impurities in noble and simple metals. 

The sharp d-bands of the transition metal impurities overlap 
in energy with the host conduction band and get broadened, 
forming virtual bound states. These become strong scattering 
centres. This idea explains fairly well the optical experi- 
ments on noble-metal-rich alloys like AgPd, CuPd, CuAg, CuMn, 
and AlMn, AlCu alloys. Again, the model merely has a 
rough tally with observations while it applies only to the 
systems listed, 

2.2 Effective Mediimi Theories; The CPa 

The effective medium theories, the t-matrix formulation, 
the Diagrammatic approaches, the Augmented Space method, the 
Recursion Method, the Method of Moments and the Cumulant 
Expansion method are attempts at refinements of the theory in 
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different directions. However, in a single site self con- 
sistent approximation, all these yield identical results - 
the Coherent Potential Approximation (CPi..)* The Average 
t— matrix Approximation (ATA) used by some authors is simpler 
than the CPA but has the disadvantage of being non-self- 
consistent. These various approaches and their limitations 
are described in the following after taking a fresh view of 
the general problem and the use of Green' s operator or the 
resolvent. 

2.2.1 The Model s 

The study of properties of solids is an example of a 
many body system. Firstly the Bom-Oppenheimer Approxima- 
tion helps separate the nuclear and the electronic parts 
of the wave function describing the many particle system 
consisting of nuclear and electrons in the alloy. The 
Schrodinger equation is still a many— electron wave equation 
describing the motion of electrons for a field s37-stem of 
nuclei. The Hartree-Fock scheme then yields a one-electron 
equation of motion for stationary states in the presence 
of the nuclei and an average potential due to all other 
electrons, say, 

(r) = (r) 

where, 

H = -■^v2+v(r) 
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V(r) being the total single particle potential. In case of 
crystals, V(r) is cell periodic and the use of Bloch' s theorem 
gives wave functions of the form, 

♦ .(?). (?) 
k k 

where u_^(r) is periodic in the real space and the eigen- 
k 

states of the Hamiltonian are eigen states of the momentum 
operator also. In case of substitutional binary alloys of 
the type 

V(r) = Z V. (r - R ) 
i ^ ^ 

where are position vectors corresponding to lattice sites 
and v^ is or Vg depending on whether the i-th site is 
occupied by the a type or the B-type constituent. The 
dependence of v^(r) on the configuration around the i-th 
site is to be ignored. Rather, the occupancies of each 
site are assumed to be statistically independent, ijffinity 
or short range order is, in other words, not taken into 
account. The probabilities of finding on A-type or B-type 
constituent at any site are then directly proportional to 
their concentrations. 

The potential V, and thus also the Hamiltonian, is 
not periodically symmetric, because the potential in each 
cell is not identical and depends upon the occupying consti- 
tuent. The periodic Bloch type electronic states will thus 
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get scattered from sites and acquire finite life times and 

« 

the crystal momentum k is no longer a good quantum number. 
Evidently, the eigenstates of momentum are not stationary 
since these are no longer simultaneous eigenstates of the 
Hamiltonian* in which case, it does not really make sense 
to consider the k-dependence of energy, that is, the band 
structure. Bansil et al (1975) have shown that it is 
possible to characterise an alloy spectrum in terms of 
complex energy bands, the real and imaginary parts of which 
correspond to the energies and life times respectively of 
the electron waves specified by the momentum k. However, 
except in the case of mild disorder, e.g. dilute alloys 
of very similar constituents, this particular approach, 
invoking as it does, integration over a complex S-space, is 
not very profitable or practicable. The method has been 
used for surface properties, but not very extensively in 
the case of disordered alloys, 

2.2,2 Role of Configuration Averaging ; 

Consistently obtained reproducible properties of 
disordered systems .in general are a surprise, considering 
the fact that each sample is only one of the many configu- 
rations possible. In particular, a realistic samj^e of 
.a random alloy corresponds to a configuration of a very 
large system. One would like to consider all possible 
microscopic configurations consistent with known macroscopic 
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parameters e.g. the concentrations of the constituents and 

take an average over these, weighted by the probabilities 

of their occurence. For a binary alloy, the number of 

N 

configurations is 2 where N is the number of alloy-sites 

in the sample. The number consistent with the values of ■ 

the concentrations of mixing is a fraction of this but is 

?0 

still too large because N 10^ even for a tiny sample. 

One could try to theoretically derive properties of the 
most probable configuration and propose that these proper- 
ties are highly likely. But again, the number of different 
configurations as probable as the most likely, and those 
of comparable likelihoods of occurrence are intractably 
numerous . 

Since the aim is to be able to theoretically under- 
stand and predict propeirties of a system, one must look at 
how these properties are obtained. These are necessarily 
results of experiments, carried out usually using a specific 
probe, e.g. photons, electrons, v^ch interact with the 
system, and are subsequently examined to deduce information 
about the system. In the process, the probe encoimters 
atoms in local environments, statistically numerous as well 
as varied. This is responsible for reproducible observa- 
tions oust as statistical systems in dynamical equilibrium 
can be unambiguously labelled by tonperature, pressure, etc. 
despite local fluctuations. The response of an alloy system. 
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for example, is thus that due to a congregation of these 
locally varied environments, the net result being an 
averaged response, not the response of an averaged environ- 
ment. This warning, that only experimentally observable 
physical quantities related to the response of the system 
to a. probe must be averaged, has been given by Anderson (1972), 
to explain source of controversies that arose over non-physical 
results obtained because non-physical quantities like self 
energy, the un— Bloch-like wave functions etc, were averaged. 

The Green’s function corresponding to any spectrxira like 
that of the electrons, phonons, etc. can be configurationally 
averaged since its imaginary part is directly proportional 
to the density of states, a meastirable quantity. Other 
physical quantities may be derived from it and comparisons 
with experimental results made. 

The above is an additional pointer to a valid local 
theory, besides the fact that long range order is absent 
while short range order is retained in the disordered state. 
The evidence of the latter is in the experimentally obtained 
field distribution function (Fig. l.l) indicating that the 
coordination number and bond lengths remain very nearly the 
same (Sec. 1.4). 

Local orbitals viz. the atomic site orbitals, or the 
bond orbitals in case of covalently bonded at cams and their 
interactions between near neighboiors only have been used as a 
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basis by chemists for small as well as very large molecular 
systems to explain the observed phenomena. In alloy systems 
too, affinities are known to affect essentially the near 
neighbour configurations. 

The case for a local theory is strong and is welcome 
because of its implication that considering configurations 
of local environments, i.e. small clusters may be sufficient. 

2.2.3 Resolvent or the Green* s Operator i 

By definition, the resolvent or the Green’s operator 
G corresponding to a Hamiltonian H is given by, 

G(z) = (z I - (2.1) 

where is the identity operator. It is defined for all 
values of z, the generalised complex energy variable, except 
on a domain viz. the spectrum of the Hamiltonian H. Realistic 
Hamiltonians are self adjoint and hence have only real 
eigenvalues. The spectrum should thus be entirely on the 
real axis in the complex z-plane. The Green's operator has 
poles at each spectral point, i.e. eigenvalue. -Projection on 
to the corresponding eigenstate is the residue at that point. 
The diagonal members G (z) of any representation of G(z) with 
respect to any linearly independent basis {|i^} is analytic 
everywhere in Ira z 0. Also, Im G^J^z) < 0 if Im z > 0 
and Im Ggg(z) > 0 if Ira z < 0. If the spectrum of H is bound 
from above, clearly Ggg(z) ~ l/z as z -* «» along the real axis. 
In case H is not bounded from above, we have to further assume 
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this asymptotic hehaviotirs of G„^(z). Finally 0(2) = G''‘(z). 
All these properties define a Herglotz function* G (z) is 

SS 

thus herglotz. 

From the definition [Eq. ( 2 * 1 )] of the Green's 
operator, its diagonal matrix element in sane representation 
{ls>} the Green's function is given by, 

Ggs(z) = < s I (2 I - s > 


_ 5. <s|n> <n|s> 

~ ^ 2 - E 

n n 

where {ln>} is the complete set of eigenstates of H, corres- 
ponding to the eigenvalues normalised to unity. Then 

the trace of the Green’s function matrix in this representa- 
tion, 

<sjn> <n |s > 

Si z - ^ 


Tr G (z) = Z Z 

SS _ j|^ 


= Z 
n 


z ^ E 


n 


T^diere the completeness of the bases set {js>} is xised. In 
case the Green’s function is required in the tight binding 
basis in a crystal, all diagonal elements of this matrix 
are identical and. 


g(z) 


n n 


It is clear now why ; g(z) is not defined at -the eigenvalues 



27 


of the Hamiltonian K, mentioned earlier. Since the eigen- 
values are all real, the g(z) is defined along the h real 
axis in the limit that the imaginary part of z becomes 
infinitesimally small, i.e., 


g(E + i6) = Lim Z 


n E + i6 - E„ 


= ^ - i7i 5 (E-E„) 


E being real, and represents the principal part so that 


- J Im g(E + iO"^) = 


6(B- V 

n 


= n (£) 

the density of states per site. This relation is very useful 
because knowing n(E), the complete Green' s function at any 
value z can be obtained from Kramers-Kronig relation. 


T «• Im , g(E+iO'*’) 
g(2) * - “• / rri;’ dS 


= / 




This relation is a consequence of the analyticity of g (z) 
and its asymptotic behaviour. 


2.2.4 Mean Field Approach; 


Here, an effective Hamiltonian is defined such that 
the Green' s function corresponding to it is functionally equal 
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to the configurationally averaged Green* s function of the 
model Hamiltonian imder consideration. That is, in a 
tight binding basis, corresponding to the exact Hamiltonian, 


H = Z ^ V, .T.. 

i ^ i,3 


define. 


H' = H + Z S.. (E) T 
■ i.d « 


^3 


where H^ is some crystal Hamiltonian, 


*-0 


^ E? + 2 Trf-i 

1 ^ i.j 

sfty, that corresponding to the pure host or the VGA, such 


that. 


where, 


H - Bo + ^ “i El + $13 


- 1 and 


and the self energy, 

s^^CE) . d (E) *0- - ?3ti i) a - 6„) 

having translational synmetry as in an ordered system# The 
corresponding Green’s operator then is 

G’(z) = (z I - H'r^ 


The idea is to seek those o{E) said o. ^ which yield, 


G’ > > 


( 2 - 2 ) 
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where G represents the Green’s operator corresponding to 
the exact Hamiltonian and the double parentheses indicate 
configuration averaging. 

In a single site approximation, the self energy is 
taken to be diagonal and identical for all sites. The exact 
Hamiltonian is then different from the effective Hamiltonian 
only in the diagonal part, 

H = H’ + Z ( w. - a(E)) P. 
i 

= H ' + h 

h being the difference and is random since is 
or (sg -s’ ) with probabilities C and (1 - C) respectively 
for a type random binary alloy. Defining now an 

operator T such that, 

G = G’ + G’T G' (2.3) 

Then, 

T = (I - hG’ h 

= Z i- p. 

i 1-d^G 

= Z t^^^ say 

i 

where d^ = -cr(E). Then (2.3) implies, 

« G » = G* + G' « T » G‘ 

For (2.2) to hold, 

« t^^^ » = 0 



This is the CPA equation, Esqjlicitly v/ritten, this is 


o 


eg - e -CT 


l-(e^- i -a)g(z-i -<7) 1 - (sg-e -<y)g(z“*e-«?) 


(1-C) = C 


The above is a single site approximation talcing into account 
only single site configiorations. Even though the effective 
Hamiltonian is still non-random, it is an improvement over 
the VGA since some information about disorder has gone into 
it. The self energy c(E) is obtained by solving the above 
equation self consistently at every energy. Physically 
what is sought is a characteristic crystjal at every energy, 

i 

which has the same Green 's function as that obtained after | 

averaging over exact single site configurations. No accoimt 
can be taken of off diagonal disorder and clusters. We shall 
make use of this basic picture of clusters (in this case only | 
one site) immersed in a mediiim which is self-consistently I 

I 

determined by the fact that no extra scattering occurs on i 

immersion ( equivalence shown in the next section) , when we 
generalise to our cluster — CPA. As such this idea is of 

I 

central importance. j 

1 

2.2.5 The Multiple Scattering Theory and Diagrammatic Techniquej 


Just as in the case of the Mean Field Approach, use 
is made of some periodic reference Hamiltonian H^, an effective 
energy dependent Hamiltonian H (E) and the exact Hamiltonian H, 
different from H by h. This difference h is the perturbation 
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causing scattering. It is the random part of the exact 
Hamiltonian and satisfies the requirement of the Multiple 
Scattering theory in that it can be decomposed as a sum of 
contributions from each site, i.e., 


h = S h . , h . = e . - o (E) 

ill 1 

Also, as defined before, 


T = h + h G’T 


= Z h. (I + G’T) *= £ T . say 

i ^ i *“ 

which, by iteration, gives, 


(2.5) 


£ h^ + Z h. G Z h. 
i i “1 “ j 


Z h, g'z h. g' Z h, + , 
-1 - j -a - ^ -4s 


( 2 . 6 ) 


which gives with (2.5) 


= 


where t^^ « 


t. (1 + G' Z T.) 

(1 - h^ g') h^ 

h^ + h G ‘ h^ + h^ G*h^ G'h^ 


(2.7) 


the atomic t-matrix representing complete scattering from 
the i-th site. Substituting this into (2.6) gives. 


T = Z t. -I- Z t. G ' Z t . + Z t. G ' Z t . G* Z t. + . 

i 1 i 1 0 .1 .^.1 k 




Each term in this expression for the total T-matrix represents 
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successive scattering processes that an electron may imdergo. 
Multiple scatterings at a particular site are included in the 
t^ as seen in eqn. (2.7). But the electron can get scattered 
again at a site* after atloast one intermediate scattering 
process » since in the third term It can be equal to i. 

As before, 

« » = G' + G’ « T » G' 

Now, writing, 

H' = IL + 2: (cr(E) - e)P. = IL+h, 

— -o i o n n ~o 

and imposing the physical requirement that the scattering 
matrix for the perturbation h^ is the same as the confi- 
g\iration averaged one for the random perturbation h, we obtain, 

H’ » + « T » (I + G'« T »)“^ (2.8) 

In the Averaged T-Matrix Approximation the « T >> 

is substituted into this to obtain the effective Hamilto- 
nian, while setting, 

« T » = 0 

gives the CPA equation. This is a self consistent approach 
while the ATA is not. 

Since the set of equations essentially solved (self 
consistently, of course) for the CPA are, 

« T.» = « t. » (1 + G« I « T. ») (2.9) 
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the approximation involved is the dropping of the fluctua- 
tion term, 

« t. G* S (T, - « T. ») » 

on the right hand side. The CPA equation 

« tj^ » = 0 Vi 

is solved for any one site because of the periodicit 3 '’ of 
the averaged quantities. 

While the ATA and the CPa yield correct results for 
various limiting cases, e.g. the low concentration, weak 
scattering and the atomic limits, onl 3 J' the CPA is reasonable 
in the strong scattering regime. The CPii has been shown to 
yield (Schwgirtz and Siggia, 1972) the first eight moments of 
the density of states function exactly and the higher moments 
also fairly accurately. As against this the jrxTA produces 
only the first four moments exactly, the RM the first two and 
the VCA, the first three. In the split band limit (|e^-£:g|»0), 
the subbands appear in the CP^^ and the aTh, but the shape of 
the curve given by the CPA is better because it yields more 
correct moments. Also both jiTA and CPi^ are invariant with 
respect to the interchange of the constituents- But the 
most important advantage of the CPA over the ATA, however, 
is reflected in the herglotz Green's function (Mookerjee 
1979) that it yields and ttie localisation theorem that it 
obeys (Kirkpatrick et al, 1970). According to this theorem 
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there can be no states in the energy region simultaneously 
forbidden by both the pure crystals of type and type B 

I 

constituents (Kumar, Joshi 1978). 

2.2.6 Other Approaches ? 

Besides the (i) Mean field approach and that of 

(ii) Multiple Scattering of Bloch states of perfect crystals 
by imperfections described above, there have been other 
independent attempts yielding the CPij. results. These are 

(iii) Similar to the multiple scattering' approach described 
above, which is a perturbative method using a propagator 
expansion, a locator formalism based on localised states 

as the bases set also gives the CPA results in a single 
site approximation. [Matsubara, Tayozawa 1961, Matsubara, 
Kaneyoshi 1966, Ziman 1969], (iv) DLagrammatic approaches 
make use of the perturbation expansion to draw associated 
diagrams, symbolising product terras of the perturbation 
terms and the unperturbed propagator, and in this way making 
systematic partial sums, e.g. those corresponding to all 
kinds of repeated scattering from a single site, which 
precisely yield the CPA, This approach was introduced by 
Yonezawa and Matsubara (1966), Yonezawa (19^) and Leath 
(1968,1970) using the propagator formalism, Leath (1970) 
also obtained the identical result using the locator expan- 
sion. 
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A formal equivalence of these and the multiple 
scattering approaches has been shown by Ziman (1969) and 
Ducastelle (1971). 

Attempts at generalisations of the CPA in the mean 
field multiple scattering and the equivalent diagrammatic 
approaches have been frustrated because of the enormous 
number of terms or diagrams required to be summed. Calcula- 
tion of the self energy over only near neighbours itself 
becomes numerically intractable. Electron self energy in a 
pair approximation using a functional derivative technique 
due to Kadahoff and Baym (1962) has been calculated by 
Schwartz and Siggia (1972). Identical results from the 
multiple scattering technique have been obtained by Cyrot- 
Lackman and DucastelleO-971). Outside the nearest neighboxor 
off diagonal self energy in any further cluster generalisation, 
it becomes virtually impossible to keep track of the terms. 

(v) A cumulaiit expansion method developed by Yonezawa and 
Matsubara (1966 a,b, 1967) to calculate the configtirationally 
averaged Green* s function, also gives the CPA result in the 
single site approximation after the correction pointed out 
later by Yonezawa (1968 a,b) is incorporated. Equivalence to 
retaining terms corresponding to repeated scattering from a 
single site has been shown. The formalism does admit genera- 
lisation to pairs, triplets, etc. but ary realistic calculation 
becomes intractable and any further approximations lead to 
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violation of herglotz property in various xmpredictable 
energy regions . 

(vi) A method of moments developed by Cyrot and co- 
workers (1970 - 1974 ) is another approach besides the cumu- 
lants expansion method mentioned above, not depending on 
band structures on Bloch’s theorem of any approximate Hamil- 
tonian. It is based on the fact that the p-th moment of the 
density of states, 

= / dE eP n(E) 

can be obtained for a tight binding model by counting the 
number of paths of length p which start and return to the 
same lattice point. The density of states can then be 
obtained from its moments. In principle all moments are 
required but , Ducastelle atad Cyrot-Lackman (1971) show how to 
get approximate p(E) knowing a finite number of moments, and 
also how to determine band edges and shapes near singulari- 
ties which this approximation smoothens. This method yields 
results closely resembling the CPA results but does not 
suggest systematic generalisations. 

(vii) Introduced by Haydock et al (1972a, 1975), the 
Recursion Method involves the generation of a new ortho- 
normal bases out of the localised orbitals and their inter- 
actions. Within this basis, the Hamiltonian is tridiagonal. 
The Green's function can then be- expressed as an infinite 
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continued fraction. The method is a purely real space 
method. It makes no appeal to the Bloch theorem or recipro- 
cal space, and as such is ideally siiited to cases T^diere 
periodicity is absent as in disordered solids, surfaces, 
solids with defects etc. Cases with more than one active 
band have also been studied by this method. Basically the 
method systematically samples the local environment of sites 
to yield an equivalent semi-infinite linear ' chain. 

The method is ideally s^iited for realistic calculations 
with the help of present generation computors. Although in 
this work the *Ray* Integration technique has been employed, 
the Augmented space formalism introduced by Mooker^jee (1975) 
to be described subsequently allow a direct use of the Recursion 
method in generating configuration averages Green’s functions 
also. The fortunate coupling of these two techniques provides 
us with the generalisation we have been seeking. 

2.3 The Aufigaented Space Method : 

Introduced by Mocker jee (1973a, b) it is an entirely 
different way of viewing the problem, but yields the single 
site CPA very naturally. This approach is powerful in that 
it can systematically be generalised to larger cluster 
apprcncimations. Essentially the Hilbert Space describing 
the true configuraticoi dependent Hamiltonian is augmented 
■fcy another Hilbert space describing the disorder, Methotis 
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applicable to ordered lattices are then used to obtain 
physical quantities. The self energy concept is not 
explicitly present in this formalism but later Mookerjee (1975b) 
furnished the equivalence by considering an electron inter- 
acting with a fermion field and derived an .expression for 
self energy identical v/ith the CPA result. 

Since the method is rigorous and tractable even when 
configurations over finite clusters are considered, this 
seems to be one of the most powerful avenues open for gene- 
ralisations. Also because of the fact that in this formalism 
one always obtains Green’ s functions that are herglotz it is 
not only attractive because of its elegance but is clearly 
superior to all other methods examined above. The method 
and its explicit form in the Cluster-CPA (CCPA) for binary 
alloys is described in Chapter III. Limitations of the CPA 
described in the following section, over come in the self- 
consistent CCPA are also discussed thereafter. 

2.4 Limitations of the CPA s 

As a single site approximation, the CPA is the best 
possibility and has been extensively used for almost all 
disordered systems of concern. It is an approximation all 
the same and yields, for example, incorrect band widths. Of 
course, no off-diagonal disorder involving two or more sites 
can be consistently taken into account in any one site 
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approximation. In binary alloys, while the centre of gravity 
of the bands is determined by the orbital energy (diagonal 
element in tight binding bases - the site orbitals) and the 
spread by interactions between neighbours, the actiaal shape 
of the density of states function is determined by relative 
values of the various interactions between near, not only the 
nearest, neighbours. These numbers can be quite different 
for different substances, even though their atomic potential 
functions are similar. 

Again as mentioned earlier, short range order is 
important in cases of binary alloys where affinities between 
the atoms of the constituents cause particular kinds of 
clusters to occur preferentially. Even in the absence of this 
effect of chemical clustering that is not reflected in the CP.'i, 
the effect of statistical clustering is also absent. This is 
believed to be important in the impurity bands i*e. when 
concentration of one component is low and 5 , the difference 
in orbital energies of the two components, is large. In fact, 
although impurity bands are produced, most of the structure 
expected in these is lost in the CPA, Machine calculations, 
however, do reproduce this structure (Payton and Visscher 1967, 
Dean 1972, Alben et al 1975), The shape of the impurity band, 
greatly affects quantities like conductivity of such impurity 
bands, the electronic specific heat, the paramagnetic open 
susceptibility, the cohesive energy (Van der Rest et al 1975a) 
etc. 
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Band tails as given by the CPA fall of much faster 
than are expected by experiments. These are understood 
(Lifshitz 1964) to be due to very large sized clusters of 
atoms. Band edges are particularl 3 ?’ important for electronic 
properties of amorphous semiconductors. 

The CPA density of states thus, in general, consists 
of roimded peaks and fast falling smooth tents. Further, 
no information regarding the extended or localised nature 
of the electronic states, and therefore also the mobility 
edge can be obtained. This information is needed for correctly 
interpreting data obtained of the electronic spectra from 
electrical and optical experiments. 

Actually, a large number of experiments concentrated 

•Ni 

alloys have pointed towards the importance of the local^envi- 
ronment on the magnetic as well as electronic properties of 
these alloys. 

The cluster CPA formalism developed in Chapter III 
accounts for statistical clusters along with off-diagonal dis- 
order over the clusters. The chemical bonding view taken 
in the application to the III-V ternary semiconducting alloys 
emjdiasises the importance of a local environment. 



CHAPTER III 


THE CLUSTER COHERENT POTENTIAL APPROXBIATION : FOI^ALISM AND 
APPLICATION TO III-V TERNARY ALLOYS 

3.1 T he Augmented StHce Method ; 

The augmented space method and the use of graphical 
methods to determine Green's function to various approximations 
is described in this section. 

3.1.1 General Theory t 

Let {e^} he the set of random numbers describing a confi- 
giaration and let P({e^}) describe the probability density 
associated with the configurations. Assuming that the various 
e^ are statistically independent, 

P({e^}) = ^ Pi (Si) 

where PiCe^^) is the probability density of the individual 
variables. Correlated distributions can also be considered. 

Siort ranged order problem has been considered by Kaplan and 
Gray (1976, 1977). 

The formalism hinges on the observation that the proba- 
bility density has all the properties of a density of states 
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corresponding to a Green’ s function of a .jiermitian operator, 

OO 

PfCe^) > 0 and / Pj^(e^) de^^ = 1 

^ OO 

so that it is possible to express this as the imaginarjr part 
of a Herglotz function which could in turn be the resolvent 
corresponding to some self-adjoint operator in some Hilbert 
space 0^, i.e., 

Pi(e.) = Lim < 4 > (^*2) 

^ ^ ^ e^e^+10* ° 1 o 

where jf^> is a member of the basis set {|f^>} spanning 0^. 

Given a p(e), if one could find a convergent continued 
fraction of the kind, 


1 



then a representation of the operator N is a tridiagonal matrix 
with a^^, a 2 , etc. along its diagonal and b^^, b 2 , ••• etc. 

along the off-diagonal i.e., 

0 0 0 . . , 

b2 0 0 . . . 

a^ b^ 0 • • . 
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in some bases Of course, we have to restrict ourselves 

to probability distributions all of whose moments are finite. 
However, almost all physical problems conform to this. The 
notable exception is the Lorenzian. The procediare involved 
in Haydock's (1972) method of Recursion for determining matrix 
elements of Green’s operator is the inverse of this. 

The question of the uniqueness of M is irrelevant for 
our purpose. Any operator M that yields the correct p(e) 
suffices. 


For the binary alloy case with diagonal disorder only, 
^c®l-c' ® represents the diagonal element of the Hamil- 
tonian in a tight binding (or 'atomic orbital) bases, and the 
p(e) is a two-peaked delta function at e^ and e^ of strengths 
c and (l-c) respectively. Setting the zero of energy at e^ 
and letting ©q “ ®B “ 

p(e) - c6(e) + (l-c) 6 (e-e^) 

And since, 

p(e) ^ Lim g(z), 

^ z-e+iO+ 

therefore, 

g(z) =■-§-+ (3.5) 




Il 
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O O 

were = (l-c) e^, - ce^t - c(l-c) e^ , 

tation of the matrix M is thus the 2x2 matrix 3 


A represen- 


(1-c) e^ 
i c Cl-c) e. 


Vc(l-c) e. 


ce. 


(3.4) 


o o 

eigenvalues of which are 0 and e^. Its eigenstates |f2> and 
lf2> describe occupancies corresponding to these eigenvalues. 
The g(z) corresponds to <fQ { jf^>, where ifQ> = 

+ (1-c)^/^ lf„> (See Eq. (3.3)). 


The ’Disorder Field’ is now introduced. To each variable 
e^ we have associated a 2-dimensional vector space 0^ and a 
representation of a self-adjoint operator corresponding to 
the probability density of that variable in some bases {ifjj^>}. 
If the set of eigenvectors say of were taken as the 

bases, then each in 0^ would correspond to one particular 

value taken by the variable e^. In general, anj'" element in 0^ 

is a linear combination of the {|h^^>}. 0^ thus describes all 

1 

possible states or configurations of the variable e^. 

The product space * = 0 . can now be constructed 

i ^ 

which contains in its description, all possible states of the 
set {e^} i.e., all configurations of the disordered system. The 
bases in the product space are, 

|f> = { jfiH}®{|fi2>}(x)... 

and if g^(z) « (zl - 


(3.5) 



then, 


P([e^}) 

= --~Im<f lG|f> 

It 0 ’ 0 

where. 


H 

O 

V 

II 


and, 

SlCz)®g2(z)®gj(z) (x)... 

G(Z) a 


$ can now be appropriately called the disorder field. All 
information about various statistical configurations is 
contained in this, while that of probability of these is 
described by the operator defined on this field. 

Configuration averaging procedure can now be described. 
Consider first, a function f (e) of one variable only. Its 
average is, 

oo 

fa/ f(e) p(e) de 

•pdcao 

= / f(e) ^ Im Lim g^^l) (z) } de 

where, 

g®(z) = <fg|{zl - 

so that, 

p(e) = - i-Im g^^(e+iO'*')» -ilm <f^|G^^^(e+iO‘^) |f_> 

Assuming now that f(z) as a function of a complex variable z 
has no singularities on that part of the real axis which forms a 
branch cut of the function g^^(z), or in its nei^bourhood, 
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f = - 2^i ^ g^^-^(z) dz 

The contour, over which the integration is carried out, is 
taken around the branch cut of g(z) along the real axis, not 
including any singularities of fCz), Since M is a hermitian 
operator on 0, we can write, 

g^~^ (z) = / (z-h)"^ p(h) dh 

—oo 

Hence, 

f = - ^ ^ f (z) dz / (z-h)"^ p(h) dh 

= - P(h) '^h ^ f(z) (z-h)"^ dz 

«aOO 

= “ ^ <fQ|G^^ (h+iO'^){fQ> dh / f(z) (z-h)"^dz 

mJOO 

= <fol[/° {Im G® (h+iO-*-)} dh f(h)]|f^> (3.7) 

..joo 

= <fol£®IV 

where f (M) is the same functional of H as f is of h, 
f(M) = / f(h) dP(h) where P(h) is the spectral projection 
operator of M. The average has thus been expressed as a repre- 
sentation of a suitably constructed operator. The generalisation 
to a function of several variables is straightforward with one 
for each variable e^. In case f(e) is itself a matrix 
element of the kind, say, <r^ {F(e) |rj>, state vectors 

I 
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H 


where, 


^ £i®%, + V t 
1 



(3.11) 


^ = I® I® ••• ... (3.12) 

The expanded Hamiltonian H contains complete information about 
the system and is defined on the expanded Hilbert space!?®)^ *, 
Description of the dynamical behaviour is in the subspace , 
while that of the statistical behaviour is in * , the dis- 
order space . The configuration averaged Green’s function is 
a particular representation in this expanded space ¥ • The problem 
has thus been reduced to the determinationof the resolvent of 

rsj 

the Hamiltonian H in the expanded space ’p , ensuring, in the 
process, the necessary herglotzicity of iiie resolvent so obtained. 
Notice that upto the final stage of configuration averaging 
no approximations are Involved. After averaging it 'is much 
simpler to introduce physically valid approximations. 


3.1.2 Graphical Methods ; 

'Etie graphical method (Haydock 1972) for determining 
various fmctionals of the Hamiltonian by renormalisation is 
summarised below. Its use in the application of the cluster 
CPA formalism for a general alloy system is made later in 
Secs. 3-1.3 and 3.1.4. 

Let {|i>} be any ccaaplete linearly independent bases 
in which the Hamiltonian is described. This may be the tight 
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binding bases^ the bond orbitals, the Wannier functions, the 
Bloch functions etc. The overlap matrix ^ 

S„ = <1|^> 

is a unit matrix in case the bases is orthonormal. Matrix 
elements of an operator A are given by, sa 3 ^, 

= <iiAjG> 

A graph is associated with every bases {|i>}. To every element ; 
li> of the bases set { |i>} is associated a vertex v^, and I 

to each distinct pair of elements (|i>, |d>) a link or bond | 
1. . which may be directional (i.e, 1.^ ma 3 ^ be distinct from i 

10 I 

1..)' or otherwise. To any operator A there corresponds a ! 

graph consisting of vertices atnd those links for v/hich A^^ ^ 0, | 
In general, the graph does not look like the lattice. It does 
so in case the chosen bases is a tight binding one, e.g. the 
Wannier functions and an operator whose off-diagonal matrix 
elements in this bases are finite only for near neighbours. 

Suppose we want to invert an operator A, Then the contri- 
bution of each vertex v^^ is defined as k(vj|^) = that of 

a link 1 . ^ as k(l. = A. ., 

Ij J-J -LJ 

A path of length N is defined as a series of ordered 

vertices and links = {v. I 4 ^ v. 3-4 4 v. ... }. 

^ o ^0^1 1 12 2 N 

Contribution of a path is defined as kCPjr) = k(vj^ ) it k(l^ ) 

n n np n p 

A self avoiding path of length N, = {v^ *^^0^1^ ^^1 *** ^ 
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is a path such that v. v . ¥ n m, with the only exception 

n 

that V. may be the same as v. in which case it is a self- 
o ' 

avoiding polygonal path. 


The renormalised contribution of a vertex v. in a path 


is defined as* 


"N 


K (v. ) = Z Z K (O 

N=0 Q'j, ^ 


(3,13) 


n 


where 0^ are all polygonal self-avoiding paths from v. and 

^n 


back, on the graph from which the vertices v^ 
been removed. 


have 


"n-1 


The following theorem then follows (Haydock 1972, Mookerjee 

1978) . 

Theorem ; 

OQ 

(A“^). . « Z Z K(CLr) where are all self-avoiding 

~ N*0 Qjj ^ ^ 

paths from v^ to v^. 

The proof is lengthy and a thorough discussion exists in 
the above two references. However, if A"*^ is the resolvent 
G s (zl - then the theorem simply reduces to the Watson- 

Greenwood renormalised perturbation expansion for the Green's 
functions G^^* 

The main purpose of the graphical representation is to 
simply visualise the extremely cimibersome algebraip expressicais 
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involved in the henotmalised perturbation expressions, parti- 
cularly when we shall be involved in generating approximations 
like the cluster coherent potential approximation. 


3 • 1 * 3 1-CPA. by the Graphical Method in the Augmented Spac e: 

Once the Hamiltonian H is defined on the extended space 
* it can be written in a site-field representation 

jnf> as, 

^,mf’ “ ^f'Snm ^ ^ff’ (3.15) 

where » is the operator in « associated with ■ttie probability 
distribution and V^, that on)^ * Ihe disorder field states are 
specified by the Set of states {If^>) occupied in all individual 
fields the exact form depending on the probability distri- 
bution p(e)* In general, one chooses the bases in 0^^ such that 
the representation of is tridiagonal so that p(e) is given 
by a continued fraction. For a bimodal density (Eq* 3*3) the 

for each h is given by Eq* (3*4) and its graph is a single- 
link chain like 


o- 


where a-i 


(i-c) Bq, ag 




a; 


b * 


®o‘ 


Finally the 


configuration averaged Green's flmCtion i& the resolvent correS* 

P' / / LIBRARY 

ponding to the Hamiltonian H (Eq. (3414). , ^ 



^ce. Afo. 
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Until now the expressions are exact. Approximations are 
made while employing the graphical technique? the exact graph, 
on which the calculation irom Eq, (3.13) is formidable, is 
replaced by a simpler graph on which the path counting is 
simpler and tractable. 

Choosing for A in Eq. (3.13) as (zl. - H), the configura- 
tion averaged Green's functions ^^^(z) can be determined by 
considering self avoiding non intersecting paths in the complete 
augmented space * . ¥alks in the augmented space imply 

the following. An electron at a site labelled by n and the field 
state jf> = lf^> Olf?><^... jf^> ... can be induced by H either 
to make spatial hops to one of the near neighbours of n with the 
matrix element V, keeping the field state the same, or it can 
remain on the same spatial site, while the field at the site n 
changes according to tie fields at all other spatial sites 
remaining unchanged. For |f> « jf^>(^ jf^> , 

G^(z) = <nf!(zl - H)'^ |mf> (3.16) 

The problem is reduced, therefore, to one similar to that in the 
ordered system at the expense of a largerspace in which the 
matrix element is to be evaluated. The superiority of this 
method over all other approaches attempted earlier is in the 
fact that approximations retaining herglotzicity of the confi- 
guration averaged Green's functions are easily made. This has 
been virtually impossible beyond the 1— CPA in all other 
approaches attempted so far. 
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Starting from a vertex labelled (o f), the first steps 
possible (Fig. 3.1) are to the near neighbours 1^^ in the direct 
space (shown by single lines) or to a different configuration 
state of the site o labelled by (o f^) by a disorder field 
hop (shown by double lines) 



Fig. 3. Is First steps of ‘walks’ corresponding to H 
in the aiagmented space. 


If Rq( 2) is the contribution of all self-avoiding non- 
intersecting paths from the vertex ’ o ’ and back in the space , 
then, the resolvent in an ordered system would be, 


P__(2) = 

°° z - %(z) 


(3.17) 


It includes contributions I ^ii^ where i are various 
neighbours of * o ’ and corresponds to the resolvent calcu- 

lated from a subgraph in which the vertex 'o’ has been removed, 
besides contributions of closed self-avoiding loops in3C . In 
the extended space, the configuration averaged Green’s function 
S”(z) can be similarly written as, 


^(z) = 1/Ez - R(z) - T(z)3 


(3*18) 
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where R(z) includes contributions entirely in the spatial part 
of the complete augmented Hilbert space ? , as well as those 
from non— intersecting paths vdiich include field hops but not 
as parts of closed loops. The latter are contained in T(z). 

The simplifying approximation is now made: Closed paths 
involving both the spatial and field hops are delinked. As a 
result T(z) = 0, Hence the graph consists of essentially that 
in the ordered system together with an extra field hop v/ith the 
off-diagonal element of in the tridiagonal representation 
(b = fcCr:^ e^ in case of binary alloys) as the link function. 
It may be noticed that then, starting from a vertex *nf* and 
hopping to a vertex 'nf^* one faces a subgraph with the vertex 
•nf« missing, which is exactly similar to the original one. 

That is, 


fi(of) _ r - F 

^of^, of^ “ ''of, of ~ ^oo 

Delinking of closed paths involving 
for this to hold. 


= G(z) (3.19) 

mixed hops is essential 


The above delinking procedure amounts of the 1-CPA 
(Bishop AMooker^ee, 1974). This is shown below. 


■pa* 

Let G'‘^( 2 ) be the resolvent corresponding to the graph 
in the augmented space after delinking in the manner described 
above, and let, 



6 


d»i+l 
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be the tridiagonal representation of the operator in each 
sub space 0^^ of $ . Then, 


G^'(z) 

= 1/[z-R^( 2 - Wj^%^(z)) - ¥j^^Gj^(2)] 



= ^00 - V 

(3.20) 

where, 



G^(z) 

= ( 3 ^^^ ( 2 )] 

(3.21) 


In the conventional 1-CPA formalism a self energy Z(z) is defined 
by, 


0=^(2) = l/[z - E{z) - (2 - E)] 

= Pqo [z - -2:(z)3 (5.22) 

CP 

where G (z) is the 1-CPA Green's function. The self energy 
Z(z) is determined self consistently frcan the 1-CPA equation, 


y p(x) dx = 0 

-«> l-G(z) (x-E(z)) 


(3.23) 


so that, 

G®(z) = f (3.24) 

where, 

A(z) = (3.25) 

Now, just as in Eg, (3.7), G^^(z) fran Eq. (3.24) can be written 
as, 

G^^(z) = <fQ|[A(z) I - M3*^|fQ>. 
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Graphically, therefore, 

G®(z) = l/[A(z) - (.W 

where, 

Gi''(z) = 1/[A(2) - (z)] 

But.A(z) from Eqs. (3.22) and (3.25) above is, 

A(z) = z - [z - S(z)] 

Hence G (z) from Eqs. (3.26) and (3.27) is the same as 
G^ (z) from Eqs. (3.20) and (3.21). The formal equivalence 
of the two-approximations has also been shov/n by Mookerjee 
(1974) using multiple scattering diagrams. 

3.1.4 n-CPA in the Augmented Space Method ? 

The essential approximation involved in the 1-CPA des- 
cribed in the last section is to ignore closed non-intersecting 
paths involving both spatial and the disorder field hops in the 
augmented space. Generalisation to an n-cluster CPA involves 
including onl 3 ?' those closed non”intersecting paths whtch involve 
disorder field hops over the n sites of the cluster. A self 
energy renormalising the interaction between each pair of sites 
in this cluster can then be defined. The corresponding n x n 
matrix S will have less than n‘^ independent components. These, 
can be determined self consistently by an iterative procedure 
starting with a diagonal form corresponding to the 1-CPA. 



Of course, "the method, will be useful onl 3 '' if tliis procedure 
converges. Mookerjee et.al. (1980) have shown that for the dlamc 
lattice with off diagonal disorder, rapid convergence is 
obtained, the band width increases and cluster features in 
the minority regions appear. 

The shortest closed path involves spatial and disorder 
field hops over two sites and consists of eight vertices as 
shown in Fig, (3.2). 



Fig. 3«2: Shortest self avoiding closed path in the j 
full augmented space involving both spatial ! 
and disorder field hops. 

3.1.5 Cluster Babedding s 

I 

Once the self consistent medium with respect to a 1 

I 

cluster of n sites is determined, i.e. for i,C} s C^, I 

(z) = (z - E^^) (3.26) 

where g|^^(z) is the matrix element of the Green's operator 
for the self consistent n-cluster meditan, Grij(z) is that | 

corresponding to some chosen crystal Hamiltonian chosen as a | 
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s'taT'tlzig point and is the matrix of self energies deter- 
mined, exact n— site clusters can now he embedded in the 
following manner. 


Let be the self consistent Hamiltonian Z Z. .T. . 

i,j 

corresponding to the n— cluster CPA medium obtained. We wish to 
consider now a model Hamiltonian which is exact over this 
cluster while the rest of the system consists of H^“^, i,e., 




Z Hi“ 

H^^^ + 
"id 


(n) . 


H' 


i» deQ 


exact 

id 


n 


Z 

i»d 


i,deC. 


/„exact p(n)N 
'‘"id ~ id ^ 


say 


n 


(3.27) 


The difference Hamiltonian SH^"^ is finite only over the clus- 
ter C^, The Green’ s function G corresponding to the If is, 


f = (2l 

■= (zl - H^"^ - 

= (I - (3.28) 

equivalent to the usual Dyson equation. Each of H^"\ 

T and ^ are n x n square matrices. The inversion of(r-G^^^6H^^^ ) 
is facilitated by the fact that 6H^"^ has non-zero elements only 
in the cluster subspace, so that the inversion reduces to that 
of a finite matrix. 


By considering clusters of successively larger sizes, 
^details of the environment implanted can be gradually increased, 
since includes exact de-tails of the interactions over the 
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cluster v^ile it sees the rest of the medium effectively • 
Choosing a central site and its set of nearest neighbours, the 
density of states at the central site obtained from the corres- 
ponding number on the diagonal of ^ is clearly closer to that 
of the true Hamiltonian than the VGA or the 1-CPA. 

Each configuration possible over a cluster considered 
can thus be dealt with to determine the densitj*" of states via 
the diagonal Green's function at the central site and eventually 
configurationally averaged Green's function determined. Parti- 
cular features of the density of states thus obtained, e.g., 
peaks, can be identified as arising due to specific clusters. 
Further, some guesses about affinities may also be possible. 

3*2 Application t o III-V Semiconducting Alloys : 

The above formalism for the n-CPA is now applied to the 
specific cases of binary substitutional alloys of III-V semi- 
conductors GaAs, InAs, GaSb, InSb, GaP and InP, These form 
zinc-blende structures. The cations (Ga, In) and the anions 
(As, Sb, P) occupy an fee sublattice each, relatively displaced 
by a fourth of the body diagdlial:" of the smallest cubic unit 
cell. Each anion has around it as its nearest neighbours, four 
cations tetrahedrally arranged and vice versa. Fig. (3»3) 
shows the connectivity of a general zinc blende structure. 

Pairs of these semiconductors with a ccanmon cation or anion 
mix in all concentrations. The sublattice of the commofi ion 
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is 'reg\ilar> or » ordered' while that of the other ion is 
occupied randomly. 



Fig* 3*35 Zinc blende crystal structure, 

3.2,1 B^ ondjng and Antibonding States as Valence and Conduction 

Bands 8 

Harrison et al (1973» 1974) have shown that bonding and 
other physical properties can be described well in the Bond 
Orbital Model (BOM) in such covalently bonded systems. The 
basis states are taken to be centred on bonds connecting nearest 
anion-cation pairs. The basis functions themselves are cons- 
tfTJicted from combinations of the sp^ hybridised orbitals of 
adjacent atoms. The bonding and antibonding states resulting 
from these broaden into the valence and conduction bands res- 
pectively when these ions are immersed into the lattice of the 
solid. 

In the BCM, the bonding functions alone are taken as the 
bases. In case the bonding and antibonding states are coupled 
significantly, the bonding bases is incomplete to describe even 
the valence band alone, as shown below. 



£v 2c "the projection operators corresponding 
to the valence band (bonding states) and the conduction band 
(antibonding states) respectively and let the Hamiltonian H 
be partitioned into subspaces of these. 


where, 


H 


mmm 


^ = 2v 8- £v > 

= Sv H Pp , 

h'^'^ = Sc 8 £v 

if = Sc8Sc 


The Green's function determined in the sub space of valence 

VC CV 

orbitals only is approximate unless H and H , the inter- 
actions between the valence and the conduction states subspaces 
are negligible, rendering the total Hamiltonian nearly block 
diagonal. 

Using the theory of partitioned matrices (Bellman, I960) 

the exact expression for the Green' s operator corresponding to 

the valence band in terms of the partitioned submatrices of H 

is derived in the following, 

P 

Let denote the inverse of an operator M in the 

sub space y, then we are seeking, 
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“ £v ^ Ey 

= ( 21 -©"-’' 

- ^ £<=> H=V)" ^ 

— P 

in terins of the partitioned juatrices, where = (z^ - H^) 

VC 

If H is small, we may to that order neglect the second term 
in the above expression. This would be true if the valence 
and conduction bands are well separated and there is little 
hybridisation. For the present computation, we confine ourselves 
to the bonding states only, ignore hybridisation and reserve the 
use of the above relation for a better calculation of the 
valence band for a future program. 


3.2.2 The BCM Bases Set ? 

The basis set to be considered are thus the chemical 
bonding orbitals. To fix the notation (same as that used by 
Chen and Sher, 1978) let us use one of the fee sublattices, 
say that of the anions as the reference lattice on which each 
lattice point is represented by a lattice vector j. The four 
bond orbitals surrounding it may then be represented by 
|'Ja> , a = 1,2, 3, 4 centred at points specified by j where 

= (1, 1, 1) a, t 2 = ^3 “ 

s= (l, -1, -l)a, where a = l/8th of s, the edge of the 
primitive cubic unit cell, as shown in Fig. (3*4), The bond 
orbitals themselves are, of course, composed of symmetric linear 





nii-N ANfi 4. 'fJsrH 



O 
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constant of the 
fee lattice 




Fig. 3.4s Connectivities and coordinates in a flattened 
diagram for a zinc blende structure • 
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combinations of the hybridised atomic orbitals ja> and |c> 
centred at the anion and cation atomic sites respectively, i*e. 


l'3a> = [la> + |c>] ^ = |ac>, say 

ac 

where N is the normalisation constant ’^2(1+3)’, S being the 
ac 

orbital overlap <a|c> of neighbouring anion-cation pair, and 


la> = I i V ^^a^ 

|c> = I js^> + 1 p^> 


(3.29) 


where |s_> and }p_> are the s and p orbitals on the anions 
a a 

while ls_> and |p > are those on the cations. The p-orbitals 
are taken to be directed along the bond in the above expressions. 


The Bloch basis functions can now be constructed from 
the bonding orbitals in the usual manner, 


l^a> 


1 , 
-=- £ e 

Vn j 


3a> 


(3.30) 


where N is the total number of anions as cations in the crystal 
and Ic the wave vector within the Brilluoin Zone* 


3*2.3 The BOyi Hamiltonian and the Eigenvalue Spectrum . 

The secular equation can, therefore, be written as, 

detlEI - iC I . 0 (3-31) 

dS 


This in the k-space becomes 
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det I E(S) I - (k) j = 0 

where, 

h““{5) = Z ^ exp [iS.(3 + - V)] 

d Od “ “ 

obtained by using the transformation of bases defined in^ 5.30)« 
One of the site suffixes on the Hamiltonian in the real i 
has been fixed as the origin without any loss of generality 
represents, for example a 4 x 4 matrix representing the 

00 

energies and interactions between the four bonds originating from’ 
the anion site at the origin. In a similar manner, assigning 
non-zero lattice vector values to j, interactions between other 
bonds can be seen to be involved, 

Chen and Sher (1978) have shown that retaining inter- 
actions upto second nearest neighbour bonds is enough to obtain 
a band structure exhibiting all required features, and that 
including the third nearest neighbour interactions, alters the 
band structure insignificantly. Ignoring ^therefore, interactions 
of further than second nearest parallel and non-parallel bonds, 
the distinct matrix elements of the Hamiltonian are e.g. of 

the kind H^ , h 3£ , H^ , H^ where - 2 (t„ - t, ) and 

00 Or^ 0r2 0r2 x 2 1 

r^ = 2( z^r ^ 2 .) T^oth symmetry translations of the fee sublattices. 
These may be renamed H^, Hp Hg, H^ for explicitness, H^(H^) 
represents the interactions between first nearest neighbour 
bonds meeting at an anion (cation) site and Hg (H^) 
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the interactions between second nearest neighbour parallel 
(non-parallel) bonds. Together with the diagonal bond energy 
say Hq, these form independent parameters describing t^ie 
Hamiltonian, and correspond to (yj +y J), (vf “ ^ 2 ' "^2 

and D respectively in the work of Chen and Sher (1978), the 
k-dependent 4x4 Hamiltonian matrix imder this approximation 
directly from Eq. ( 3.32) is, therefore, 

“ HQ+2H|[cos(2x+2y) + cos( 2 x+ 22 ) + cos( 22 + 2 y)] 

= HQ+2H|[cos(2x+2y) + cos(2y-2z) + cos(2x-2z)] 

H^^ = H|^+ 2 H 2 [cos( 2 x- 2 y) + cos( 2 y- 22 ) + cos(2x+2z)] 
H^ * HQ+ 2 H 2 [cos( 2 x- 2 y) + cos(2y+2z) + cos(2x-2z)3 

H^^ = 2{H^cos(x+y) +2H®cos(x-y)cos2z}-2iH®sin(x+y) 

H^^ = 2{H®cos(x+z) + 2H^cos(x-z)cos2y}-2iH^sin(x+z) 

H^^ ss 2{H®cos(y+z) + 2H^cos(y-z)cos2x}-2iH2sin(y+z) 

H^^ = 21H®cos(y-2) + 2H^cos(y+z)cos2x}+2iH^sin(y-2) 

H^^ = 2{H®cos(x-z) + 2H^cos(x+z)cos2y}+2iH^sin(x-z) 

H^^ ss 2{H®cos(x-y) + 2H^cos(x+y)cos2z}+2iI^sin(x-y) 

where H® » (H^ + H^)/2 and h| » (H^ - H^)/2 and x, y, z 

symbolise ^ k^a, ^ k^a, ^ k^^a respectively. 
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3.2.4 The Parameters ^ 

The values of the parameters in Table 3.1 have 
been obtsiined by conveniently fitting energ 3 ' separations 
at specific symmetry points on the Brilluoin Zone, r , X,L,etc. 


Table 3.1s Distinct matrix elements of the 

Hamiltonian used as parameters for 
the III-V semiconductors between 
nearest neighbour bonds and second 
nearest neighboiir bonds. 


Mate- 

Parameter 

rial 


H? 

*4 



GaAs 

-1.5187 

- 0.4250 

0.3406 

- 0.1094 

InAs 

-1.3937 

- 0.5125 

0.3844 

- 0.0656 

GaSb 

-1.2687 

- 0.3750 

0.3906 

- 0.0594 

InSb 

-1.3500 

- 0.3875 

0.2937 

- 0.0563 

GaP 

-1.3625 

- 0.3250 

0.3812 

- 0.1188 

InP 

-1.3125 

-0.4125 

0.2625 

- 0.0375 


The parametrization approach employed is empirical. 

Also fitting with experimental data from different sources 
gives different results. Often variation in the values of a 
parameter obtained by fitting to different data for a given 
material is more than that over different materials, Chen and 
Sher eventually impose an energy dependent Lorentzian broaden- 
ing and vary the peirameters in order that the broadened 
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peaks 'thus obtained best resemble the experimental spectra 
obtained by Ley et al (1973). 

The above procedure is somewhat arbitrary and its 
validity is also unclear. Despite this, the values of the 
parameters have been used for a preliminary cluster-CPA cal- 
culation. Results of these are presented in this chapter. 
Dissatisfaction with the above parametrization approach 
motivated an ab initio cbemi'coX pseudopotential calculation 
for the required parameters, starting with Herman-Skillman 
tables of self consistent Hartree-Fock atomic orbitals and 
potentials. The procedure for this calculation, comparison 
with the empirically obtained values and the results of the 
cluster-CPA calculation constitute the follov^ing chapter. 

The numbers shown in Table 3.1 are fitted with the experimental 
data by Chen and Sher (1978), However, if we fit to each expe- 
rimental data, the numbers are different. For example , for 
GaAs, the difference in the parameters for the same material 
is more than that between different materials (Table 3»2). 


3.2.5 Matrix Elements oi the Green* s Operator in the VGA and 
1-CPA s 

The alloys under consideration are of the kind 
Ga^In^^^^s as a common cation alloy, and. GaAs^P^^^^ as an 
example of a common anion alloy. These alloy in all propor- 
tions and therefore x may anything between 0 and !• Each of 




Table 3*35 Bond lengths, bond energies and work 

functions for the III-V Semiconductors 
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the six materials GaAs, InAs, GaSb, InSb, GaP and form 
the zinc— blende structure of nearly the same bond length 
(Table 3»3)» We shall assume, therefore, that the interac- 
tion parameters corresponding to thfe pure materials are 
the same as those in the alloy. 

Top of the valence band as the zero of the energy 
scale used by some authors is incorrect. The vacuum state 
must set the zero of the energy scale since that is the 
common state. Different -wrork functions for the various 
materials (Table 3.3) also supports this viev/. 

The VGA parameters are given by the average of 
those of the constituents, weighted by their respective 
concentrations. For dilute alloys, the parameters corres- 
ponding to the host are also good enough as a zeroeth app- 
roximation for the later CPA calculations, but in these 
cases, the virtual crystal is very nearly the host anyv/ay. 

Matrix elements of the Green* s operator for the virtual 
crystal corresponding to the alloy are sought. An efficient 
method due to Chen (1977) called the ’Ray* Integration Technique 
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is employed for the purpose. Essentially a systematic k-space 
integration over the Brilluoin Zone (BZ), the general method 
is described in Sec, A, 1 of Appendix A, The 3-dimensional 
integration is reduced to far simpler one-dimensional inte- 
grations over thin tetrahedra into which the irreducible part 
of the BZ is divided. The method is fast and attractive for 
a machine calculation. 

-Fig. 3. 5 shows the density of states for pure In^ as 
obtained by the ‘Ray’ integration method. Band structure, 
i.e, energy eigenvalues along scane standard symmetry direc- 
tions of an fee lattice is also shown for the same material. 

In the form presented by Chen (1977) the method is 
applicable only to the calculation of functions whose inte- 
grands in the k-integral bear the symmetry of the BZ, The 
diagonal Green’ s. function can, therefore, be obtained by 
direct application of the method. 

In Sec. A-2 of Appendix A, a generalisation of the method 
applicable also for functions, whose integrands do not reflect 
the symmetries of the BZ, is presented. The off-diagonal 
matrix elements of the Green’s operator are in this class of 
functions and have been obtained by this generalized *Ray* 
integration method. 



0 . 




Fig, 3.5* Band stmcture of InAs and its density of 
states "by *Ray’ integration technique 
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Ideally, a self-consistent calculation should be 
carried out. This has been stressed time and again by 

authors v;ho deal with the coherent potential approximations. 

A C P 

However, since the variation in the parameters , 

is rather small as compared to that in the bond energy 
Hq over the various materials (Tables 3.1 and 3.3), cluster 
embedding is proposed to be carried out in the relatively 
easily obtained self consistent 1-CPA medium rather than the 
self consistent n-cluster-CPA medium. The use of a 1-CPA 
medium instead of the self consistent one, as is done by 
Kaplan and Gray (1976) leads to a narrower band width and 
distortion of the density of states near the band edges. 

The self consistent calculation on a single s-state on the 
diamond and bcc lattices has been recently done (Kumar, 

Mookerjee and Srivastava 1980» 1981 , Mookerjee and Srivastava, 
1981). However, the full self consistent calculation on the 
diamond lattice with foursp^ hybridised orbitals per site or 
the four bonds per pair of sites is still beyond the scope 
of the present work and will have to be left for later research. 
For the present cases under consideration, the corrections 
are expected to be small, but for other alloys with large off- 
diagonal disorders the effect of self consistency may be crucial. 
Under the above approximation, the self energy is a scalar, 
determined by solving self consistently eq. (2.4) reproduced 


here. 
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The VGA diagonal Green's function is substituted as a zeroeth 
approximation and the quadratic equation in a thus obtained 
is solved. The correct solution is chosen by checking that 
a) the chosen solution goes to zero as 6 , the difference 
the bond energies of the constituent materials, goes to zer 
and b) that herglotz property is maintained. The Green’s 
function is then calculated at E - o and substituted back 
into the above equation. The process is iterated until self 
consistent solution is obtained for a fixed value of energy. 
The energy parameter is spanned across the valence band. 

Convergence is fast at all points except for those 
around sharp peaks in the VGA density of states. The real part 
of the self energj?" crosses a zero around these regions, and 
therefore the value of the self energj?" is taken to be the 
average of that at neighbouring points. 

Fig. 3.6 shows the real and imaginary parts of the self 
energy, and the VGA and the 1-CPA diagonal Green's functions 
for ^s'x^l-x^® ^ 0*3, 0.5s 0,7, 0,9. 

3.2.6 Embedded Clusters : 

Eq, (3.28) is now directly applied. A choice of the 
clusters considered is, of course, to be made at this stage. 





Pig. 3*6* (1) -Ira (11) -Im (lit) -Ira i: and (Iv) do t 

tor the specific alloys Qa^lnj|^_^jjAs with x « O.l, 0,3# 
0*5* 0*7 end 0*9 obtained from the fitte»l parameters. 
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For a cluster of n sites, G®, G^^\ jjC n) unit matrix 

I involved in the equation eire n x n matrices. 

Figs. 3.7 (s) to show various distinct configura- 
tions of each cluster and the corresponding probabilities 
of occurrences are drawn in the following. Sites are indi- 
cated by small circles and bonds by lines connecting these 
sites. The alloyed ion is indicated by coloured circles. 

Sites are labelled for the specific alloys ^a^x^l-x:^* 
which the cation sites are alloyed and for which grajdiical 
results are presented subsequently. Exact clusters for any 
other allojr can easily be obtained from these diagrams. 

The n X n difference Hamiltonian 
_ j^exact __ j^Cn) 

where is the 1-CPA Hamiltonian as an approximation to 

the self-consistent n-CPA medium Hamiltonian. A particular 
cluster may consist completely of only one of the alloying 
constituents. In that case H^^ simply consists of inter- 
actions between the i-th and the o-th bonds as in the 
imalloyed state. For example, for the alloy Ga^^In^^^^s, 
one of the configurations of any cluster under consideration 
will consist of only Ga and As ions i.e. only Ga-As bonds. 
Depending cjrt whether i,;3 are identical, nearest neighbour bonds) 
parallel second nearest neighbour bonds, or whatever, the 
corresponding interaction parameter of pure Ga-As crystal 

03CSLC*t# 

would be substituted in Hyj • 



(a) 


(c) 



llg. 5*7 U)f O) ^ ©ontinuad. 

I 

t 
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The first term is vanishing because of the orthogonality 
of two different hybridised orbitals at the same site. The 
middle two terms would be identical in the unalloyed case due 
to the geometrical symmetry of the location oi the orbitals 
involved. Hence this part of the expression is equivalent 
to the average of these corresponding to the pure cases. The 
last term is small compared with the middle two since it 
involves more distant orbitals. The approximation, therefore, 

involves the neglect of the variation in this term over the ; 

1 

constituent materials, which is small anyv/ay. 

i 

Various clusters considered are described in the [ 

I 

following t I 

I 

(a) For a one-bond cluster, only two distinct configurations j 

are possible. The matrices become scalars and, therefore, the | 

M 

cluster Green’s function G cori-esponding to each of the confi- 
gurations and shown in Fig, 3.7 (a) are calculated 

directly using Eq, (3.28) and the weighted average of these 
is taken. Here invblves only the difference of the 

specific bond energies of the GaAs and the InAs bonds from 
the 1-CPA value. The 1-CPA Green' s function is reproduced as 
is expected, 

(b) A cluster of four bonds meeting at an alloyed site is now 
looked at. In case of Ga 2 j.In^_^As, the central site is a cation 
site and there are again only two distinct configurations possib: 


k 
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since only one alloyed site is involved in the cluster. These 
are indicated in Fig, 3.7 (b) with the four bonds numbered as 
shown, the 4x4 matrices are 


Ho 


Gq gJ g= of 

4 Hq 4 4 

. = 

4 °o 4 4 

T-T ijC 

“"0 % 


4 4 4 

ttC ttC ttC tj 

^0 


4 4 4 °o 

^ ^ 





C 

where Hq and are the exact bond energy and the nearest 
neighbouring bond interaction parameter for bonds meeting at a 
cation site. Gq and are, in a similar notation, the diago- 
nal and the appropriate off-diagonal matrix elements of the 
Green’ s operator in the bond basis, in the 1-CPA. G^^^ is the 
same for the two configurations and shown in Fig. 3.7(b) 

while in H , appropriate numbers for the specific bonds and 
interactions are plugged in. Again, Eq. (3-28) is used to deter- 
mine the diagonal Green's functions corresponding to each cluster 
and an average taken weighted by the probabilities x and (1-x) 
for the two clusters, to yield the configxjrationally averaged 

, M 

Green's function. The diagonal elements of «G » so obtained 
are identical due to symmetry and correspond to the S’ sought - 

(c) A four bond cluster in -vdiich the central site is the 
unalloyed site (the anion site in case of Ga^In^_^As, for exam- 
ple) is considered next. Its specific configurations and 
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As a starting reference, a series of Ga In, As den- 

X 

sity of states per bond scaled by ir in the VGA and the CPA 
are shewn in Fig. 3*6 together with the real and imaginary 
parts (with the sign changed) of the diagonal CPA self energy 
for various values of the concentration x indicated. The 
curves reveal an s— like band separated by a gap from a p— like 
band. The latter also is apparently made up of two separated 
peaks. These correspond to the second nearest neighbour 
parallel and non-parallel bond interactions. This is substan- 
tiated by the observation that the integral under the second 
peak is twice that under the first one. 

As expected, the CPA roiinds off the sharp peaks present 
in the VGA, lowers the 6 -function like spikes and has small 
at low concentrations of either constituent, increasing to a 
maximum at the 50-50 alloy. 

Results of the density of states per bond scaled by it 
(that is, the imaginary parts of the diagonal Green's functions 
along the real axis with the sign reversed) for each configura- 
tion, and that averaged over these, for each of the clusters 
Cl, C^, C^t, Cy and Cg described above (and shown pictorially 
in 'Pig. 3.7(a) to (e)) are shown graphically in Fig. 3.8 for 
the alloy OaQ^^ In^ g As and in Figs. 3.11 for the alloy 
Gao 5 ^^0 5 pairtial density of states corresponding 

to different clusters are important in that they indicate the 
way the effective medium spreads out the cluster eigen-energies 





Fig# 3*3(a),(b) and (c) continued. 




Pig* 3*8 1 Density of electron states per bond in units of 

(-Db 0) for the alloy ^,9^* correapondiag to 

each configuration (Fig* 5.7(a) to (e))of the olusters 

(a) (b) (c) (d) Cy and (e) Cg and those • 

averaged over oonfigurations of eaeh olusteri tisiai 

the fitted parttseters* 
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Fig. 3.10 (a) and (t) continued. 




P4g* 3*10 (o) and (d) oontiauad. 
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?lg. 3*10t 


(a) -.lm« 0 »« , (to) -Ia«a»« , (c) -IiK< 0 »« , 

* 

(d) »X»« 0 »Q^ and (a) -la«Q» 0 ^ for 0 «o,l ^** 0*9 ^ 
using tha fittad pannatars* 
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Fig, 5,X3 (a) and (ti) aontlimad. 
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Obq^jj Ai 



Pig* 3*131 (a) -Hb «0»q , (b) -Ib«0»q , (c) «lli«0»g ^ , 

1 A A 

(d) -!■ aiid (a) -fli«0»g^ for Q«Oi5 ^0*3 

using tho fittad iwraBatars* 
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when the latter are immersed in the former. The str'jcture 
of those partial densities also indicate the cluster origin 
of the structure appearing in the averaged densities of states. 
For example, the enhancement of the peak at the top of the 
hand in the 7-bond cluster CPA for GaQ ^ Iuq ^ As as compared 
to the 1— CPA and 4*— CPA can be attributed to the clusters 
Cy^^^ and Cy^^^ whose peaked structure at that energy contri- 
butes to this cluster property (See Fig. 3.6 and 3.10), The 
cluster effect is evident pronouncedly in each of the clusters 
consisting of either constituent of the alloyed ion only. 

Partial densities of states are also of interest in the calcu- 
lation of the optical properties, particularly, the structure 
near the valence band top. In the absence of other calcula- 
tions to compare with, it is necessary to generate these 
partial densities to account for the difference between the 
CPA and the cluster-CPA’s. 

Fig, 3.9 shows the VGA and CPA density of states on a 
larger scale for a more detailed examination. Fig. 3»10 (a) 
and (b) show that embedding single bond clusters is insufficient 
in case we want to describe these ternary alloys, vSiile embedd- 
ing the four bond cluster with the alloyed cation site at the 
central site, reproduces results resembling the 1-CPA very 
closely. C^ is, therefore, the smallest unit embedded that 
is self consistent. Again, comparison of Fig. 3*10{b) and (c) 
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shows that in case of a four bond cluster, the one* with the 
central site alloyed or disordered gives more information 
than with the unalloyed site at the center. 

Fig. 3.10(c) and (d) compare the 4' -CPA and the 7-CPA. 
The 7-CPA involves embedding a central bond with its six 
nearest neighbours, and this can take into account variations 
in and and some of the second nei^bour bond inter- 
action. The significant clustering effects arise at the band 
top. This can be understood as due to the fact that only the 
states at the top of the band are relatively extended enough 
to sample the environment in detail. Lower lying states are 
localised and so do not sample the environment in detail. 

Fig. 3 *10(8) shows the result for a 6 bond ring 
cluster Cg embedded in the medium. The results are compar- 
able to the 7-CPA. The second near neighbour interactions are 
sampled but not completely. The cluster variations which 
give rise to the enhancement of the upper peak are also 
present here. 

The Fig. 3.11 presents the partial densities of states 
for Gan K Ihn K As. The cluster effect is less prominent now. 
This is to be expected as in model calculations (Kumar, 
Mookerjee and Srivastava, 1981), The partial densities also 
indicate that there are no highly probable peaked structures 
in these clusters. The CPA smooths out the VGA (Fig. 3*12) and 



S9 


■fch6 cl'us'bsr CPA s further introduce sniootheningj lov/ering 
peaks and broadening out the density of states. There is 
again a perceptible filling up of states at the band top, 
if we compare Pig. 3.13(c) and (d) as we go from the 4’ -C.FA 
to the 7-CPA. The ring cluster (Fig. 3.13(e)) is as before, 
not very different from the 7-CPA. The clustering effect , 
though smaller than in case of x = 0.1, can be clearly seen 
tc increase gradually as the size of tfte cluster is increased, 

3*3 The Impurity Bands t 

The appearance of impurity bands in a model calcula- 
tion is a triumph for the 1-CPA, The n-CPA cluster embedding 
in a model calculation for a 90 - 10 alloy for 6 = 1.0 
has been done by Mookerjee and Srivastava (1980) for the 
s-band. The graph (Fig. 3.14) clearly shows the extra lumps 
appearing. The cluster embedding in it is reported to broaden 
the band perceptibly. 

3.4 Conclusions s 

The n-cluster CPA formalism developed above is a 
refinement of the 1-CPA preserving essential physical 
properties of the configurationally averaged Green’s function, 
gradual change in which to be expected as the size of the 
cluster embedded is increased. 

The application of the above formalism to the valence 
bands of III-V semi-conducting ternary alloys, indeed, shows 
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"the above • The clustier effects in these alloys is not 
dramatic# This is due to the fact that the variations in 
the diagonal energies are small compared with the* widening 
faclsor , where is "the syimne'bric corabinaiion of the 

nearest neighbour bond interactions, and Nevertheless, 

the effects are not negligible specially in dilute alloys. 

That the environment is not self consistent, is a drawback 
and may mask some further clustering effects, specially 
because of the off diagonal disorder. The self consistent 
problem could not be handled because of difficulties both 
in the formulation [specially in deciding the self-consistency 
condition in the cluster to the interior and sm'face atoms 
(Butler, ( 1973 ) ] and in implementation of some of the pro- 
posed schemes, viiere the equations involved too many 
variables and seemed intractable. Yet even with the non-self 
consistent medium some features of clustering do show up. 

This in itself is the achievement of this work* Similar 
embedding in a non-selfconsistent media have been tried on 
model systems (Kaplan, and Gray, 1976, 1977, Kaplan et al 1980) 
and have shown that considerable amount of cluster structure 
does show up. However, it cannot be emphasized enough that 
eventually the self consistent solution should be obtained 
in future . 
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■fch© six III— V SGifliconductors dealt with in the last chapter. 
The values of the parameters are compared with those obtained 
by parametric fitting. 

Thereafter, in Sec. 4.4 the VGA and the 1-CPA medium 
calculations are carried out. Finally exact clusters are 
embedded arid Green’ s functions corresponding to each specific 
cluster and that obtained after a configuration averaging 
over these clusters is done, are presented graphically for 

4,2 The General Chemical Pseudopotential Scheme r 

The fact that simple metals have energy bands that 
are free electron like is surprising because the conduction 
states are required to be orthogonal to the core orbitals 
and must therefore be quite different from plane waves. The 
reason for this is that the Schrodinger equation for the 
conduction electrons can be transformed into a pseudo- 
Schrodinger equation which has the same conduction band 
eigenvalues as the original Hamiltonian but for which the 
wave functions are slowly varying in space and are nearly 
plane waves. Only for S— vectors near Brilluoin Zone boun- 
daries may these be combinations of degenerate plane waves. 
Clearly, there is no unique transformation that alone makes 
this possible. One such pseudopotential equation is the 
Austin et al form (1962) 
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T}0^> + (V - E |0 ><0 JV) |0^> = E 1?! > 

k c ^ ^ ^ S' 

The summation is over the core orbitals. Part of the poten- 
tial V is cancelled by the core terms on each atom so that 

the net potential is weak. The eigenfunctions are thus 

k 

little distorted from those of the kinetic energy operator 
T, i.e. the plane waves. 

Similarly Anderson (1969) introduced a pseudopotential 
such that the pseudov/avefunctions look localised and atomic 
like. The formulation Justified Huckel* s scheme for parametri- 
zation of molecular interaction in a fashion clearer than the 
usual procedure of making a unitary transformation to an 
orthogonal basis. 

The aim of localised orbital theories is to present 
the exact band molecular orbitals say, as linear combi- 
nations of orbitals |0 > localised at the atomic sites. In 

UEr 

case the distortions due to neighbouring orbitals is small, 
this provides a chemical building block picture (Boys I960, 
Boys, Foster I960, Gilbert 1964) in which the wave functions 
of large molecular systems are built up from chemically 
invariant parts, which need not be calculated for each new 
system. Also, correlation effects can be conveniently studied 
in this basis. These ideas can be applied to crystals consi- 
dering these to be very large molecular systems. 
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The Wannier functions as the basis set, are not 
suited for the above picture since the orthogonality cons- 
traint between those associated with different sites necessa- 
rily produces oscillations, more dispersed than the atomic 
functions* Also, the tails of the functions are rather 
sensitive to the environment. It is now clear that better 
sets of localised functions can be defined by employing the 
pseudopotential techniques (Adams 1961, 1962, 1971a, 1971b, 
1971c, 1974, Gilbert 196h, Anderson 1968, 1969). 

Consider a system made up of N atoms with a given self- 
consistent one-electron Hamiltonian H. In the isolated state 
these atoms have well defined energy levels and eigen- 

uu 

functions |0°> due to the free atom Hamiltonian such that 

Ha|02> = e°|0°> (4.1) 

When brought together to form a molecule or a lattice, these 
atomic orbitals overlap, and in the latter case, foim bands, 
each containing exactly N states? 

We are looking for a useful set of local orbitals 
linear combinations of which yield the s. Clearly these 

span the same band space, hence, 

IV = piv 

where P is the pro^Jection operator on the band i.e. 


(4.3) 
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N 

P = E I 

i=l ^ ^ 

Using "tlie properties of commutation with the Hamiltonian and 
idempotency of the projection operator, i.e, 

PH = HP and P^ ?= P 

Eq. (4.3) can be rewritten as 

H|0^> - PHP|0^> = 0 (4.4) 

In order to define 10^> s uniquely, another condition must be 
incorporated. The most useful j0 > s will be those that 

wu 

resemble s most closely. To derive the j0„>, v/e divide 

the total system into subsystems v/hich may be atoms, ions, 
bonds etc, and constrain t0„> to satisfy H , the Hamiltonian 
associated with one such sub-system as well as possible within 
the permitted function subspace, i.e. 

Within a single band, N nontrivial functions |0„ji> will satisfy 
this Eq. (4.5), but only the lov/est one will be 

locali^sed and will correspond to the atomic like orbitals of 
Ha (Bullett, SSP, Vol. 35). Adding Eqs.(-4.4) and (4.5), we 
obtain, the Adam 's equation, 

H|0^> - P(H-H^)IV = 

or HJ0^> + {U^ - HJ 3 P}| 0 „> = e^|0„> (*.6) 

where 5 H-Ha, is the perturbation to the cell Hamiltonian 
due to the rest of the system. 



This is the determining eqmtion for j?i^> ejid shows 
why the environment may often have little effect on the 
original localised, orbitals of Eq, (4.1) since the pertur- 
bation itself is reduced by projecting out all ot it which 
acts within the band. The residual screened interaction 
- PUqP can most often be treated by standard perturbation 
theory by an iterative procedure to obtain self consistent 
orbitals 10„> starting from 10°>, However, Anderson noints 
out in his work (1969), that whenever the scheme is applied, 
no attempt to solve for the exact eigenfunctions is required, 
that these can be found or exist in principle is all that is 
necessary. 

It may be noticed that the localised orbitals corres- 
ponding to different cells are determined by different Hamil- 
tonians in general, and there is no orthogonality relation 
between them. 

The pseudopotential equation introduced by Anderson is 
obtained by simplifj^ing the projection operator P to 
involving projections only on the local states and 

U = Z V, 

^ b;^a ^ 

where is the perturbation due to cell labelled by b. 

Then using the properties PP^ = Pl^ and P|0jj^> - 10£j^* Eq. (4,6) 
becomes. 
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l.e. H|0^> = V ^ I IV <"'•7) 

Anderson applied this to the case of the Ti-orbitals of 
benzene and showed that the cancellation is so complete that 
the self —consistent solutions are ver^?’ close to the unpertur- 
bed atomic orbitals. 

The molecular orbitals are given by linearly combining 
these self consistently obtained orbitals, i.e. 

IV “ ('‘•e)- 

The usual secular equation determining the eigenvalues would 

then be 

det iH - ESl = 0 (4.9) 

where S is the overlap matrix, not a unit matrix since the 
basis is not orthogonal. A far more efficient and convenient 
method which takes full advantage of the form of the Eq.(4.7) 
would be to rewrite it as, 

H|0„> . Z |0p> (4.10) 

where and Dp^ = <0p \0^> for h-4 a, and substitute 

this in the Schrodinger equation to give 

HiV = ZC„iH|0„> 

- ^ I l’'p> 

I V IV 

a 



Now linear independence of the orbitals on different sites 
alone implies that the coefficient of everjr j0 > must vanish. 

Cc 

Hence the secular equation becomes simple? 

det I = 0 (4.11) 

It may be noticed that orthogonality of the bases 

is not demanded, nor does the overlap matrix appear in the 

esqpression. Further, the off-diagonal matrix elements D . 

( \ 

contain only the localised operator which is relatively 
insensitive to modification by the neighbouring cells. 

The price, of course, is paid by the non-hermiticity , 
in general, of the pseudo-Hamiltonian matrix obtained since 
<0p and <0gj j0p> may not be alike unless each 

cell is the same and is placed in the same environment. 

In fact, from Eq. (4.10) 

D = S“^ H 

and therefore, D will be Hermitian if and only if the Hermitian 
matrices S and H commute. This will be so when every orbital 
environment is identical so that these are diagonalisable 
simultaneously. Yet, no imaginary components should arise 
in the eigenvalues since Eq. (4.11) is trivially related to 
the Hermitian form in Eq. (4.9). 

Premultiplying Eq. (4.6) by <0^j we obtain. 
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In C3.S6 "bh© unp6r'turb6<i orbitals are not very different 
from the self-consistently obtained orbitals [0 > from Ea.(4.6) 
the former can be substituted for the latter in the above 
Eq. (4,12) to yield 


= 




(4.13) 


giving the orbital energy in the bonded state in terms of 
that corresponding to the unperturbed state and a first order 
correction. The correction consists of the first order pertur- 
bative term lessened by the orbital overlap times the inter- 
active term. 


In case the bases chosen consists of the bond orbitals, 

an effective bond Hamiltonian H „ associated to a bond orbital 

ac 

|ac> needs to be defined such that 

Hg^^|ac> = E^|ac> (4,14) 

where is the bond energy and the bond orbital |ac:> consists 
of the site orbitals ja> and {c> on the bonding atomic sites i,e,, 

lac> « jj^(|a> + |c>) (4.15) 

ac 

where N is the normalisation constant » 2(1+S), S = <a,c>, 

ac 

the overlap. The orbitals ja> and jc>, in case of diamond 
and zinc blends sti*uctures are the sp^ hybridised orbitals* Now 


H 


ac 


» T + V, 


ac 


writing 
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4.3 Application to III-IV Semiconductors t 

Just as before* interaction terms upto second nearest 
neighbouring bonds are retained, this time of the pseudo- 
Hamiltonian. The pseudopotential interaction terms to be 
calculated are those, between adjacent bonds with a cation 
or an anion site in common €uid the second nearest neighbour 
parallel and the non parallel bonds. A direct volume inte- 
gration is carried out in a 'Wigner-Seltz * cell defined 
around a bond. This consists of a three dimensional space 
such that every point in it is closer to the centre of this 
bond than the centre of any other. 

Fixing the origin of the coordinate system on the centre 
of the bond, and the z-axis along the bond in one direction, 
coordinates of the centres of the neighbouring bonds and those 
of the sites involved are determined, A grid is assigned to 
the spherical coordinates (r, 6, 0) such that the iiiree dimen- 
sions of the curvilinear volume elements so obtained are 
comparable. These are then considered in successively larger 
spherical shells. The distances of these from the afomic 
sites concerned are used to determine the values of the atcanic 
potentials and the radial parts of the s- and the p-orbitals 
from Hermann-Skillman tables of these for the cations Ga, In 
and the anions As, Sb, P. Since these are tabulated at discrete 
points, linear extrapolation is carried out between near points. 
The ' cos © » factor for the p-orbital wave functions corresponding 
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•fco 6&ch sit© is dsierminsd. ty a. scalar product of appropriate 
unit vectors. Volume elements are then considered one by one. 
In case the considered one is closer to the centre of the bond 
and origin of the polar coordinate system than any other 
bond centre, the contributions in it to the integrands for the 
parameters is accumulated. Besides the matrix elements between 
the nearest neighbour bonds meeting at an anion and a cation, 
which are symmetric, both kinds of each of those between second 
nearest neighbouring parallel and non-parallel bonds are calcu- 
lated. In a self-evident notation, these are 
D^NP(A)^ corresponding to <l|Dj2>,<llD}3>, 

<1|D|^>» <L|Dp>,<l|Dp> and <llDl7> respectively with the bonds num- 
bered as shown in Fig. 4.1. The procedure enstires integration 
only over the Wigner-Seitz cell around the bond without going 
into the details of its geometry. 


The bond potentials for Ga-As along the bond and along 
lines parallel to it in the plane containing a neighbouring 
bond at either end are shown in Fig. 4,2. 


The values of the six parameters Dt® = 


I,^NP(A) ^ 

computed above for the six III-V semiconducting materials 
considered in the last chapter acre tabulated in Table 4.1, 
along with the corresponding parametrically fitted values 


for comparison. 
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Pig* 4.1 1 The neighbourhood of a bond for 
irr-V Semiconductors. 
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arlong bond 

displaood by Q«1 of bond length 
displaced by 0,2 of bond length 


Fig* 4«2i Bond pot 'Atials for GaAa along the bond and 

along 1' iies parallel to the bcNod at a diatanoe 
0,1 and 0,2 of the bond lengfch| In the plane 
oontainlag a ndlghbourlng atom at either end* 



Table 4.1s Ccanparison of Parameters Obtained by C3iemical Bond 
Pseudopotential with those Obtained by Fitting. 
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The paraiiie"ter "that essentially gives the ’spread* to 
the bands is the symmetric combination + D^^)/2. 

The values obtained from chemical bond pseudopotentials are 
smaller than the corresponding fitted parametric values. 
Otherwise there is no other systematic trend. The parameter 
^ (D^^-D-|^^) , which measures the asymmetry between the 
cation and anion composing the III-V semiconductor, is obtained 
of different signs for different materials. This feature is 
to be expected but is absent in the values tested by Chen and 
Sher, 

If interactions upto first nearest neighbouring bonds 
only were retained, the sign of vrould not alter the 
eigenvalue spectrum (See Sec. 3.2.$ but when interactions 
between mere distant bonds are taken into account, the sign 
of becomes important. Since the pseudo Hamiltonian is 
non-hermitian, the matrix elements between second nearest 
neighbouring parallel and non-parallel bonds each are of the 
two kinds. For the purpose of the calculation, the symmetric 
combinations and are used to define an approxi- 

mate hermitian Hamiltonian, and the eigenvalue spectrum 
evaluated therefrom. 

Fig. 4,3 shows the band structures and the density of 
states of GaAs and InAs as obtained by using the fitted 
parameters and those obtained from <diemical bond pseudopotential 


parameters 



fa) GaAs 


(fc) In As 
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Fig. 4,3s Band structures and density of states per bond 
in units of for (a) QaAs and (b) InAs fron 

(i) Chemical pseudopotential parameters and 

(ii) the fitted parameters* 
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The principal difference is in the lower s— like band. 

In the LPK and XTX sectors the branch is convex towards T in 
the fitted parameters, but concave in the chemical pseudo- 
potential method, giving rise to (a) a much diminished gap 
between the s-like and p-like states, and (b) a contraction 
of the total band width. The upper branches are not very 
dissimilar. 

It must be emj^iasized, as in Bullett (1975) 
that the chemical pseudopotential method, at least in the 
relatively crude, non-self consistent calculation is quantita- 
tively inferior to the fitted parameters. But it is a starting 
point of an ab initio derivation of the energy spectrum, 
taking into account the non-orthogonality of the basis. This 
is an important factor, in case there are a variety of diffe- 
rent experimental results and we are not sure v/hich we should 
select for our fitting, 

4,4 CCPA Calculation Using the Chemical Bond Pseudopotential 
Parameters i 

The cluster-CPA calculation has been made again, now 
using the chemical bond pseudopotential parameters. The VGA 
and the 1-CPA density of states per bond times tc and -Im Z 
and Re E along the real Z axis for Ga^^ ^l-x x =s 0,1, 

0,3, 0.5, 0.7, 0.9 are shown in Fig. 4.4. The one bond 
cluster C^, the four bond clusters and C^* with one and 
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(d) 


(e) 







T 


Range 

mbsoiasat 

-20 aV to 0 «V 
ordinate I 
0 to 3 



—1 

g» 4,5t Density of eleotron states per bond in units of n 

(-la Q) for the alloy 0*0* X ^**0.9 ^ oorresponding to 
eaoh configuration (Fig. 3*7(a) to (e)) of the clusters 
(a) (b) C4 (c) C^t (d) C7 and (e) Cg and those 

averaged over configurations of each cluster* usins 
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^ 0.9 ^ 



Range 

abflolaaat -20 eV to 0 eV 
ordinate t 0 to 3 



fig* 4. 6 1 (1) -la 0^^ and (11) -2n 0^**^ for 

Q»0.1 ^®0,9 waing eheniaal bond 
paeudopotentlal paraaetera. 




Pig* 4«r (•) «nd (b) eontlnueii. 
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^* 0,1 ^ 0.9 



1 



(d)' 


Fig* 4,7 (c) «adi (d) eontinutd. 
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®* 0*1 



Pig. 4.7 1 (a) -Tin«0»p , (b) -1««0»- , (c) «I»i«a»« , 

^4 ''4 

(d) and (a) -1 *«G»q^ for lno.9^* * j 

using ohaaleal bead paaudopotantlal jauraBatars. 



?l8« 4*8 (a), (b) and (o) continutd. 












Fip?* '4*3* Density of electron states per bond in units of % 

(-Im G) for the alloy G«o«5 ^“o.5 corresponding to ! 
each configuration (Fig* 3ii7(a) to (e)) of the clusters | 
(a) (b) (c) C^t (d) Cy and (e) Cg and those [ 
averaged over conflgtff’Stions of each cluster, usliiK ! 



rig. 4«X0 (c) and (d) continuad* 
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I 





Fig. 4.101 («) >Xh«0»^ 

1 

(d) 

^7 


, (h) ^imi<Q»Q f (tt) -XM<a»Q » 

4 4 

•ad (•) -1««0»Q ivr Aa 

Q 


using ehenioal liond patudoptfttatiaX ^ariattars. 
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4.5 Conclusions t 

In conclusion it must be said that despite the simpli 
city of the definition of the bond pseudopotentials, the 
"Wigner trick" and the free atom self-consistent orbital wav« 
fionctions and potential functions used, the results are highi 
satisfactory. The main achievement being that arbitrariness 
involved in a parametric fitting with experimental results, 
which considerably differ among themselves, has been by-pass4 
The only parameter used in the application of the chemical 
pseudopotential scheme above that has its source in experi- 
mental results is the bond energy with respect to the vacuian 
state. The part of this which is the work fmction, is a 
parameter obtained very easily by a simple experimental set-uj 
for the photoelectric effect. Hence its value is reliable. 
The "bond energy" is a parameter that enters every local 
chemical bonding behavioxar of atoms and is, therefore, well 

known. 



CHAPTER V 


LOCALISATION 


5.1 Introduction : 

The question of the nature of electron states in a 
disordered system is now tackled. Unlike the extended Bloch- 
type wave functions in crystals or ordered systems in general 
Anderson (1958) first demonstrated the presence of non-diffu- 
sive ie localised electronic states. For a simple tight- 
binding model Hamiltonian with diagonal disorder and nearest- 
neighbour interactions only (now known as the Anderson 
Hamiltonian) he considered the evolution of an electron 
initially in one of the tight-binding basis states. In case 
of an oitiered system} the electron would eventually diffuse 
over the entire system and the amplitude of the final state 
at the original site would become vanishingly small. He showed 
that for some minimum degree of disorder (measured, say, by 
the width of the probability distribution of the random 
members on the diagonal of the Hamiltonian matrix) , the 
amplitude at the site of origin becomes definitely (i.e. 
with probability * 1) finite. After some controversy over 
this work of Anderson (Ziman 1969, Lloyd 1969, Brouers 1970) 
localisability became an established fact (Thouless 1970, 
Anderson 1970, Mott 1970). For various models of one- 
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dimensional systems, Ishii (1973) has shown that axr/ finite 
amount of disorder makes all the states localised. 

On making a more detailed examination, Mott (1967, 1968) 
introduced the idea of sharp mobility edges separating the 
localised and the extended states in the bands, and showed 
plausibly, that as disorder is increased, these mobility 
edges move into the bands, until these meet. He termed this 
the Anderson transition , i.e. when all states become loca- 
lised. The same picture was also : later arrived at by Cohen, 
Fritsche and Ovshinsky (1969) in their model for amorphous 
semi-conducting alloys. Now known as the Mott-CFO model, it 
is consistent with Lifshitz^s (1964) arguments concluding that 
for mild disorder, band edges tail and a finite portion of 
the band near the edges consists of localised states. 

Economou and coworkers (1970a, 1970b, 1972) define a 
"localisation function"to locate the mobility eciges and obtain 
the above movement of the mobility edge with increase of 
disorder. 

In this chapter a delocalisation criterion for the 
electronic states based on the size of the localisation 
domain' in introduced and discussed. An averaged localisation 
length is obtained. Results are presented for we^ disoixier. 
Near the mobility edges, these are soluble, while near the 
band edge, numerical computation is required. 
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Mooker^Jee (1975) has examined the localisation cri- 
terion of Abou-CSiacra et al (1973). and pointed out that the 
criterion involved quantities like <|aj|^(E)|> where a^(E)ar© 
the tight binding amplitudes of the localised wave function 
^ (E, r) at the site r^. It was also indicated that a better 
criterion should rather involve quantities of the kind 
<}a^(E){^> since it is these that are related to ftoysical 
probabilities. In the localised regime, this implies the 
study of averaged size of the dcanain over which the wave 
function is localised. This is related to the ’ localisation 
length’. Such a study would have two-fold signific^ces 
firstly, we take as our delocalisation criterion , the pro- 
babilistic divergence of the localisation domain, and thus 
determine the mobility edge E^ as a function of the disorder 
parameter. This stuc^ has already been reported by Kumar 
et al (1975) for the weak disorder limit* We shall also 
indicate the analysis in the case when disorder is larger, 
leading to Anderson transition* Secondly, in the localised 
regime, we can look at the behaviour of the averaged size of 
the localisation dcmain and the way in which it diverges as 

E - E . The size of a localised state is an Important 
c 

parameter that enters expressions for various physical 
properties, for instance, pionon— assisted hopping conduction. 
(Austin and Mott 1969)* This calculation of the averaged 
inverse size of the localisation domain will foins the bulk 
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of this present work# We shall present results for the 
weak disorder limit. 

5.2 Preliminary Formulation : 

Let us follow Anderson (1958) in conducting a hypo- 
thetical experiment. Suppose at t » 0, we place an electron 
at the site r^^ and look at the probability it is there at 
t -♦ «». This probability 


P * / dE lim |Gqq (E + iTj)|^ (5.1) 

■ n -*0+ 

Gqq(Z) is a representation of the resolvent (ZI - H) ^ of 
the Anderson Hamiltonian H « 2 ^ 

can be expressed as a Feenberg perturbation series (Feenberg 
1948), 


Gqo^z) 


oo 

2 


2 V. 


n»l 0^ 

(O, in ) 

®i i 
^2 2 


Oi, 


V . 

^ 1^1 ^ 1^2 


(O, i-j^, ■ . 


^i i 
n n 




(5.2) 


where ^ are all non-intersecting paths of length n starting 
and ending at r^. Subscripted index indicates that the 
ygj;»'tex is deleted end the Green* s function is calculated 
on the resulting subgraph. 
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V/e may now introduce the Cayley Tree Approximation 
(CTA) in tv/o vrays. Graphically, we may think of it as 
delinking all closed non-intersecting paths, i.e. ignoring 
contributions of all such paths to Green’ s function. As 
an example, in the exact lattice, there was a contribution 

^01 ^11^ ^12 ^00 ^ ^23 ^*^30 shov/n in Fig. (3*l)> 

On delinking, such a contribution is absent. 



Closed ’ path' Delinked 'path* 

Fig, 5.1'' An example of a »path» (i) in the exact 

lattice and (ii) in the Cayley Tree Approxi- 
mation, obtained by delinking, 

A study of Eqn. (5.2) shows that apart from the 
n = 1 and n = 2 terms, all the rest comes from such closed 
non-intersecting paths. Equivalently, the CTA may be 
thoxi^t of as truncating the Feenberg series after n = 2, 

The latter view has been taken by Abou Chacra et.al. (1973)* 
The first view point perhaps gives scane insight into the 
electron diffusion experiment. An electron on a Cayley 
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Suppose K is the connectivity of our Cayley tree. Then 
the self energies have the following simple relationships. 

, (Z) - 


a^(Z) = 


(K) 

Z 

keN. 

J 




IfZsE + if} (t)>0) and a ^ ~ ^ 3 * 


EjCZ) 


^j(Z) 


keN. (E-ejj-Ek)'^ + (ri+hk)' 


and 


V^(t} + 

T 


keN^ + (t5+ 


(5.4) 


Now let Tj O’*". In the localised regime, 
and hk ■* 0, except on a set of points of measure zero. If 
we replace all equalities with almost equality (almost 
= except on a set of points of measure zero), then 


E. 




» £ 1 

keN^ 

£ 

keN. (E-ek-Ek)‘ 


(5.5a) 


(5.5b) 


Note that while -«» < E^ < •», 0 < 1^ < •» because of 
herglotzicity of cf(2). Assume now that ISf s and l*s caa 
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either side of the equation have identical distributions, 
and make the substitution 


X » 


V" 


Y . 

X 


v;here, ^ = x 

1 I- > y 

which implies that f(x) = E-Sj^-x = E^, and 

g(y) = y-i = ij^. 

Now, if G(x, y) is the probability density of x, y and 
F(Ej^, Ij^) is that of Ij^, then 

•(•CO 00 

G(x, y) = / dsj^ / dE^^ / dlj^ F(E^, ly) 

— “ o 


6[Ej^-.f(x)] 6[li^-g y)3 P(ejj) 

400 

= / dej^ F(E-ej^-x; y-1) pCsj^) 


+« yZ yy2 y4 

So, G(X, Y) = / dej^ F(E-ej^ - ~ l) ^4- pCe^^) 


X 


Then (5.5a) and (5.5b) give, 

F(k, s) = [G(k, s)]^ and 


FQ(k, s) = [G(k,s)3 


K+1 


(5.6) 


that is. 


FQ(k,s) = CF(k,s)3 


(K+l)/K 
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where F(k, s), FqCJs, s) and G{k, s) are the Foiirier-Laplace 
transforms of F(Ej, 1^), Fq(E^, 1q) and G(X, Y) respectively. 
These yield the integral equation; 

r t f ■rr2 

F(k, s) = { / dx / dk p(k ) F(k , -- 2 ') 

mJCSO <mCO X 

exp [-ik’ (E-x) - 3^^ (5.7) 

where p(k) is the Fourier transform of p(e). 

We shall attempt a solution of the type, 

F(k, s) * f(k, {x^}) P(s, {Xj,}) (5.8) 

Here {x^} denotes a set of parameters describing the pro- 
bability density of E^. We note that the probability density 
of E. can be obtained directly from (5.5a). It is not 
affected by the probability density of 1^ which is why we 
have kept f(k) independent of the parameters describing 

P(l.). On the other hand, the variables 1^ and E. are 
0 u J 

dependent. So P(l^) can be affected by the probability den- 

U 

sity of Ej. This is the assumption involved in writing 
F(k, s) as in (5.8). The assumption does not imply the 
statistical independence of E, and 1.. It is an approximation 

J u 

none-the-less. Examine the particular form taken for F(k,s), 


F(k, s) 


f(k) P(s) 


(5.9) 
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for p(k) . e -^sgnCk'y k 

try f(k) = exp [-iZ^k - s(k) k] 

P(s) = exp [-as-ts*^], a>0, 0<a<l 

as a solution. It is clear from the solution that a, a, b 

are functions of Z^, However, v/e are also asvuning here 

a and b to be independent of k. Eventually we need only 

a(k =0), b (k = 0) but a(k), b(k) v/ill satisfj^ integral 

equations similar to AAT equations (Abou Chacra et.al. 1973) 

for which a full knowledge of a(k), b(k) for all k's is 

required, in order to solve for a(k = 0), b(k = 0). The 

exact import of this approximation is still mclear. However, 

it is not as drastic as assuming E., 1. to be completelj’ 

J J 

statistically independent, neither is it the severe upper 
limit approximation. 

Suppose pCS},) has a Cauchv distribution 
and try the solution, 

- 2. sgn (k) k 
f(k) = e ^ ^ 


then (7) gives for s = 0, 
iZ^k-Z^sgn(k)k 



exp [- sgn (k* )-iZj,-Z^ sgn(k' ) 
+ i(E-x)] k’ 3^ 
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. {i f dx 

^ (E-Z^-x)2 , (y->2,)' 




exp [-ikV^/x]}^‘ 


Therefore , 

+ iZ, = 


KV^ 


1 ^ ~ — it Y + 

Again, for k = 0 (5.7) gives 


P(s) = { / 


(Y+Zjj)A 


(E— Z^— x)2 + (y-»-2^^2 


? (4) 


exp (- 4)}’' 

X 


If FgCk, s) = ^0^^^ ^ q ( s ) 3 then from (5.6) 

foM = 


and 


Po(s) = [PCs)] 


(K+l)/K 


sV^> 


( 5 . 10 ) 


(5.11) 


Because of the factor exp (- '=^), the integral in 
(5.11) cannot he carried out in the usual manner. If we want 
to study the case for any value of the disorder parameter 
then P(s) must he numerically solved from (5.11) and then 
Pq(s) obtained from (5.12). 

5.4 The Weak Disorder Limit s 


If (Y + Z^) « 1, then 



1A9 


(y+Z.)/T p „ 


Neglecting the subsequent terms is equivalent to evaluating 
the integral in (5.11) by the steepest descents approximation. 
In this approximation, (5.11) reduces to 


PCs) - exp (- --SSl- ) [ ? [ 


(E-Z^)' 


(E-Z ’)2 


K 


Let, 


P(s) = e"®®Q(s) 

substituting this in (5.13), we have 

KV^ 


a = — — 


if 


(E - Z^)^ - KV^ 


Scs) = {3[— 


(E-2J2 


(5.13) 


(5.14) 


(5.15) 


If Q(x) is the Laplace inverse of Q(s), then (5.15) shows 
that it must be a stable, infinitely divisible distribution 
(Gredenko and Kolmogorov, 1954). Further, since we must have 
Q(x) = 0 for X < 0, an approximate form for Q(s) is 

exp (-s°^ Sec ^ ), 0 < a < 1. 


Substituting this in (5.15) we obtain, 

2a 2a 

(E - Z^) = KT 

therefore, 


a 



In K 



In (■— ) 


( 5 . 16 ) 
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Hence (5.12) gives. 


Pq(s) * exp [-a^s - ^ (Sec 2^- ) s“3» 
K+1 


^0 *“ lt“ ® 


The required probability density of L then is 

Po(y) = Qq 

where, 

1 '’X+iCO y. , T 

“ ScT ? T X 

X— ioo 

Oncp ^(y) is known from here, we can calculate, 


<1* - f Uq (y - ag) dy 


“ I ®o<y) 


(5.18) 


Equations (5.10) and (5.14) give the mobility edge. If a -• «, 
since PQ(y) * 0 for y < a, the domain size diverges in 

probability. This happens as, 

(Ec - Z^)2 = KV^ 

which, for y « 1 gives, 

E = 2 K V - (fKV)^/^ Y^(K^«nar et al, 1975)' 
o 

Thus E^ -* Ejj as y -* 0, correctly, as opposed to the AAT 
limit. The new criterion involving <|a^(B^j^> has the 


e(3ge over the other. 
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(i) 


Let us examine the form ^(s) = exp [- Sec s^] 

jy . T ^ T 


Since a < 1, - 


0 

Ts 


„ K+1 JX-1 r- f 'i o^/v 
= a s Qq(-; Sec 


oo 

as s - 0, since a < 1. That is, the average f y Qrv(y) dy 

0 ^ 

does not exist, as remarked before. However <1/L > exists 
and that is why v/e usxaally speak in tems of <l/L>”^ rather 
than <L>, 


(ii) If QQ(y) is the inverse Laplace transform of 

• -s“ ' 

have the same form except 

that the x-axis is scaled by a factor Sec This 

is apparent from (5.17) by changing the variable s to 
(|^j- cos s’ , It is, therefore, enough to invert 

_ga 

e and then simply scale the inverse. 

OO 

<L"^> = / 1 9^{y) dy (5.19) 

^0 y(E^ Sec + Hq + 1 

Some (^(y) for E not very near are shovm in Fig, 5.2. They 
have similar shapes, reminiscent of Pearson curves. They are 
non-zero for positive arguments and have long tails. As E "♦ E_, 
the curves peak near 1,0 and become narrow, but retain the 
tails caused by 0 < a < 1, although here the tails fall 
off faster. 

(iii) T^ftien E is very near E_, a becomes almost 1, and the 

c 

Qo(y) '^^(y-l)? or, equivalently, 
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Fig. 5.3 shov;s the averaged, inverse localisation length 
in the localised regime (E^ < E < Eg) for a very weakly dis- 
ordered system, specified by f = 0.01. The points near E^ 
are evaluated from Eq. (5.20) while those near Eg are evaluated 
from Eq. (5.18) by integration. The linear divergence of the 

localisation domain E marks the sharp transition to extended 

c 

Bloch type states in the interior of the band. 

The mobility edges act as a band edge for transport 
properties depending upon mobile electrons in extended 
states. Gaps as measured by optical and electrical experi- 
ments do not tally for this reason. Phonon assisted hopping 
conduction (Austin, Mott 1969) is possible in the energy 
range E^^ < E < Eg. This is distinguished from that due to 
extended electrons by the temperature dependence of the 
conductivity . 

5 . ^ Conclusions 8 

A quantitative estimate of the averaged localisation 
length for v/eak disorder has been achieved. The linear 
divergence as the mobility edge is approached, is obtained 
as expected. Also, the mobility edge is located numerically. 

In generalisating to strong disorders we are faced 
with a very difficult mathematical problem. TheS -function 
approximation to the Lorenzian has to be dropped and the 
non-linear integral equation cannot be reduced to a non-linear 
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algebraic one. The mathematics of the solution of such non- 
linear equations is still in an eaibro3^cnic stage. This led 
us directly into the interest in non-linear problems reported 
in Appendix C . However, that problem is a rauch simpler 
one and although it yields some insight is insufficient for 
the localisation solution. 



CHAPTER VI 


I 


CONCLUDII«J REMARKS 

The development of the theory of alloys beyond the 
1-CPA is significant, even though in the cases of III-V 
semiconducting ternarj'' alloys, to Vfhich it has been applied, ' 
the additional features obtained are not very large. But the , 

i 

path is paved for application to alloys of materials so 

widely different in mass, potentials and bond lengths that | 

1 

the VGA or even the 1-CPA are not satisfactory. The n-cluster, 

I 

CPA density of states for each specific cluster can be very , 

■ ! 

useful in determining affinities by comparing the configu- j 
rationally averaged density of states with the experimental 
one and identifying some configuration or configurations 
required in more or less than statistical properties, for a 
good agreement. 

In conclusion, we may sum up the achievements of this 

work« 

1) Formulation of a cluster extension of the coherent 
potential approximation for alloys vdiich is (a) herglotz at 
all energies, concentrations and degrees of disorder, and 
(b) tractable for realistic calculations. 



CHAPTER VI 


CONCLUDING REMARKS 

The development of the theory of alloys beyond the 
1-CPA is significant, even though in the cases of III-V 
semiconducting ternary'' alloys, to which it has been applied, 
the additional features obtained are not very large. But the, 
path is paved for application to alloys of materials so i 

widely different in mass, potentials and bond lengths that . 

1 

the VCA or even the 1-CPA are not satisfactory. The n-clust^ 
CPA density of states for each specific cluster can be very 
useful in determining affinities by comparing the configu- ; 
rationally averaged density of states with the experimental 

I 

one and identifying some configuration or configurations 
required in more or less than statistical properties, for a 
good agreement. 

In conclusion, we may sum up the achievements of this 

work,. 

1) Formulation of a cluster extension of the coherent 
potential approximation for alloys which is (a) herglotz at 
all energies, concentrations and degrees of disorder, and 
(b) tractable for realistic calciilations , 
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2) Qiemical pseudopotential determination for the electri 
nic structure of the pure materials and the alloy. 

3) Application of the above formulations to a series of 
III-V ternary alloys, 

4) Formulation of a localisation criteria based on the 
averaged localisation length, and determination of this para- 
meter in a weak disorder limit on a Cayley tree. 

! 

* 

5) Eigenvalue and eigen- functions are calculation for the 
quartic anharmonic oscillator in all regimes of anharmonicity 

(k - plane) . The accuracies of these is restricted only by ; 

1 

the machine precision. The method evolved is fast even in | 
those regimes where previous calculations have been done. 

1 

Several refinements and extensions over the calculation 
presented in this thesis suggest themselves. Some of the 
immediate ones are the following: 

1, Band structures of the III— V semiconductors show no 
dispersion along the Z and X axes of the Brillouin Zone even 
if the third nearest neighbour interactions in the Bond Orbita^ 
Model are included (Chen and Sher 1978). This may be overcome 
if the complete non-hermitian pseudopotential Hamiltonian 
matrix is solved for eigenvalues rather than the approximate 
hermitian matrix. The antisymmetric combinations of the 
somewhat different interactions between second and also the 
third nearest neighbouring parallel and non— parallel bonds 



156 


each, must be included. As pointed out in Sec, 4,2, the 
eigenvalues of the non-hermitian pseudo-Hamiltonian must 
still be real since these are the same as those corresponding 
to those of the true Hamiltonian by definition (See Eqs,4.9 
and 4.11). 

1 

2. The evaluation of the n-cluster n-CPA self consistent 
medium •which was by-passed for application to the III-V semi- 
conducting alloys must.be carried out. However, the size of 
graphs in the full axigmented space increases very rapidly. For 
a cluster of n-sites, the size of graphs to be considered are 
2n'^l_gon • Also, the number of distinct self energy functions 
corresponding to each distinct matrix element of the Hamiltonian, 
and hence the number of equations to be solved simultaneously 
for self consistency increases, making the process tedious. 

A 2-CPA self-consistent calculation has, however, been reported 
since (Kumar, Mookerjee, Srivastava 1980), 

3, The pseudopotential scheme can be used beyond Anderson’s 
approximation to include the effect of the core orbitals. Of 
course, finally the more localised pseudopotential wave functions 
and eigenvalues must be determined self consistently, although 
the use of the free atom orbitals and potentials together with 
the truncation involved in the Wigner trick is a good starting 
approximation , 
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4, Just as the bonding orbitals have been used for evalu! 
ing the Green's functions for the valence bands, the anti— 
bonding orbitals, in principle, determine the conduction banc 
However, a scheme similar to the BCM is not expected to 
reproduce the conduction band. The conduction band states 
are necessarily more extended and hence interactions between ’ 
rather distant bonds are required to be included. 

The effect of the conduction band on the valence band ; 
can, however, be studied by incorporating an approximate ; 
Green' s function corresponding to the subspace of the conduc-' 

1 

tion states, e.g, in a free electron model (See Sec. 3.2,1). | 
This is expected to alter essentially the top of the band , 
again, since these states are the nearest to the conduction 
states on an energy scale. It may be interesting to compare ; 

I 

the effect of the cluster-environment with that due to coupl- 1 

I 

ing with conduction states, since both seem to alter the top j 

i 

of the valence band, important for all transport properties. | 

5. The variation of the bond length has also been neglec- | 
ted. In the pseudopotential scheme, the interactions can be i 
re-calculated at the inter-atomic distances corresponding 

to specific alloys. In this way, the approximation made 

for Ga-As-In kind of near neighbour bond interactions can also i 

i 

be refined. This involves the calculation of somewhat more ' 
number of interaction parameters, but it is nonetheless j 


tractable 
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6. The study of cluster effects in case of vddely 
different constituent materials showing a separated impurity 
hand in the 1-CPA, can also be made in a model calculation. 

7. The calculation of averaged inverse localisation length 
for large disorder has been attempted, but the difficulties 
with the mathematical formulation are enormous. The case of 
weak disorder becomes tractable essentially because of the 

6 -function approximation to a highly peaked Lorenzian. As 
soon as this afiproximation is dropped, the already complicated 
integral equation becomes very difficult to handle. 

8. The quantum anharmonic oscillator dealt with in 
Appendix C (Banerjee, 1976, Banerjee et al, 1978) has, since 
the work was published, made the solution to aaiiy <ybhor per- 
turbative problems that are divergent but as3mptotic possible. 
The difference equation method, in vdiich scaling of the bases 
'is easily incorporated, is thus more powerful than believed 

until now. Besides the eigenvalues and eigenfunctions, all 
moments between any two levels can be obtained starting with 
only the first one as input to an appropriate moment recurssion 
relation (Banerjee 1977). At higher values of the energy level 
quantum number, the V/KB approximation can be shown to preserve 
the scale (Banerjee 1978a). The method has been extended to 
study a general inharmonic oscillator (Banerjee 1978b) , in 
which case a rescaling of the perturbation series has been 
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APPENDIX A 


GENERALISED ’RAY’ INTEGRATION TECHNIQUE AND APPLICATIONS 

The *Ray* Integration method developed by Chen (1977) 
determines essentially the diagonal Green’ s function via the 
density of states function. A generaO. spectral function 
F(E) is of the form 

F(E+) = R(E) - 1% 1(E) = f d\ —£^1 

BBZ [E'^-E(k) ] 

where E"*" means that the function is to be evaluated by includ- 
ing an infinitesimally small imaginary part in E, and E(Ic) is 
the band structure. For f(lc) = 1, 1(E) reduces to the density 
states? 

1(E) « / d\ [E-E(k)] 

BZ 

since, 

eW •4“-i^6(E) 

where E and e are real, liftien f (Ic) = 1, or more generally, when 
it is invariant under the symmetry transformations of the 
Brillouin Zone (BZ), it is clearly sufficient to integrate 
over the Irreducible part of the Brilluoin Zone (IBZ), i.e., 
that part of the BZ which, through the application of these 
transformations, reproduces the con5)lete BZ. For those spectral 
functions for which f (k) does not reflect all the symmetries of 



the BZ, as is the case with the off diagonal matrix elements 
of the Green’s operator required for the CCPA, the method as 
described by C3aen is not directly applicable. For such spec- 
tral functions, a generalisation of the method has been deve- 
loped which does not necessitate integration over the complete 
BZ, This is described in this iippendix after a brief descrip- 
tion of Chen’s procedures recapitulated for the sake of comple- 
tion. 

A. 2 Chen’s Rav Integration Method ! 

First, the IBZ is decomposed into one or more tetrahe- 
dra as triangular prisms. Consider now, one such tetrahedron 
PABC defined by the three vectors q2, q^ (Fig.A.l). Any 
point P inside it is identified by a k vector. 

1 & = aq^ + aPq2 4- ap q^ 

where a = rA'/rA, p = A’L/A’B’ and y - LP/LM, each parameter 
ranging from 0 to 1 to scan the tetrahedron completely. A'B’C* 
is a plane parallel to ABC passing through P and IM is parallel 
to B’ Cl Points corresponding to constant values of a describe 
planes parallel to ABC, Constant Y planes are all hinged at 
the edge r'A and vary from TAB to rAC. Similarly, constant § 
planes all pass throu^ r and span the tetrahedron from FBC 
to rA, cutting the face ABC at a line parallel to BC, p =0 
corresponds to the line TA. Along any line starting from 
called a ’ray’ by Chen, p and Y are constant, while a goes from 
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0 to 1, This is the advantage of the above manner of parame- 
terisation since it facilitates the use of one-dimensional 
interpolation procedures, e.g, by calculation E(a^) on a 
reasonable grid {a^} along some fixed 'ray' , interpolation 
by say, Lagrangian formula. 

E(a) = Z L. (a) E(a, ) 

can be used, where L^(a) are interpolation coefficients depend- 
ing on the order of the interpolation rule. 


The tetrahedron is, therefore, further divided into 
thin tetrahedra (TT) as shown (Fig. A. 2). Let I^CE) be the 
result of integration over the i-th TT, Then, 


1(E) = I I.(E) = Z / ds. dk f(k) 6(E-E(k)) 

i ^ i ^ 

If the TT are thin enoijgh, angular dependence can be appro- 
ximated as constant, and, 


I^(E) ~ A / dk^ kj f(k^) 6(E-E(kj^)) 


where, ar^^ 


^i 


-1 


dE(k. ) 

= AR, zf(kp i-wr\ 3 

r 1 k^=k^ 


solid angle element of the i-th TT, 
is along the axes of the i-th TT, 
root of the equation E-E(kj,) = 0 


and the summation is over various roots of the same. If k^ is 


the length of the axes and k^ = dk^, then. 



= 3V (AS^/s)/k^ 

= 3V6^/k3 


where V is the volume of the whole tetrahedron, 6^ = 

being the area at the end of the TT and S the area of ABC. 
Then, 



where = k^k^. 

The method can thus be summarised in the following: 


(1) liLvide each tetrahedron into thin enough TT's specified 
by andYj^ corresponding to their axes, i.e. line Joining 
to the median of the trian^Aar face on the base, (ii) Determine, 
E(a) at a few points n in number along the axes of 

each, including one at the tip r and one at the BZ surface, 

(iii) Fix the energy value E at which the value of the imaginary 
part of the spectral function is to be determined, (iv) Use an 
n-point interpolation formula, say Lagrange’s, to evaluate the 
values of at which the value of E, if at all, is hit b 3 r say, 
Newton-Raphson technique, in each TT, considering all of these 
one by one. (v) Determine the slope dE(a)/da at these points 
by calculating E(a) at points ^/2 around s and dividing by 
ha, choosing ha to be small enough, (vi) Evaluate I^(E) as in 



contributions in each IBZ,-like portion of the BZ is to be 
evaluated. As shovm in the following, these contributions can 
be obtained by the use of the symmetries of the BZ, without 
having to solve the eigenvalue equation again. 


Suppose a point Ic = (k , k , k ) has been located in 

Y H 

the IBZ at which the value E is hit in the n-th band, that is 
E = EjjCk), then, corresponding points in the six similar IB 
Zones (three of opposite parity) in the quadrant x > 0, y > 0, 
z > 0, are obtained by permuting the components of k in all 
possible ways, i.e. the six points are described by coordinate 

sets XCkjj, ky, ^ 2 ^* ^^x* ^y^* ^^z* ^y* ^x^* 

(ky k^, k^) and (k^, k^). Changing signs of k^, ky and k^ 

takes one to the octant described by x < 0, y < 0, z < 0 and 

permutations of these gives corresponding points in the six 

IB Zones in that octant* Similarly, changing signs of one or 

two components at a time and permuting each such set, the rest 

of the six octants are scanned. This identifies the 48 points 

in the BZ at which the value E is crossed, corresponding to the 

point (k , k , k ) located in the IBZ actually scanned. 

X y z 


Writing the point k as a 3 x 1 coliimn matrix 




k. 


y 


each of 


the above transformations can be written as a 3 x 3 matrix which 


when postmultiplied with the column matrix yield another column 
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matrix representing the corresponding point in one other 
IBZ. For example, the matrix. 


interchanges the values of x and y components of k. Similarly 
Sy 2 defined. Together with the products S^Sy^.* 

®xy^x 2 identity operator, these form a six-membered 

group under the matrix multiplication operation. That is, 
products of any number of these in any order yield a matrix 
identical with some member. The product is, in general, non- 
commutative and the group is isomorphic to the symmetric group 


0 10 
10 0 
0 0 1 


S 


z 


Further defining the operator S 


-10 0 
0 10 
0 0 1 


which 


changes the sign of the x component only, and similarly, S^y, 

S , their products between themselves and those with the six 
defined above, cover all the 48 parts of the BZ similar to IBZ. 
These 48 matrices also constitute a group as can easily be 
checked. 


The Hamiltonian from Eq. (3.32) is given by a 4 x 4 
matrix at each 5 

= 5 iC' [exp ±{lS.(? +7 (A.3) 

r rO 



A-9 


where aa* are bond orientation indices and the sian is over 
vectors connecting sites in the same fee- sublattice. It will 
now be shown that applications of any of the transformations 
described above, amounts to renaming these indices a, a* • 

Writing the column matrix corresponding to the point 
^ as K and the displacement vector R = (r + ), as ^ 

k.R = K^R 

If k'= S K , where S is one of the 3x3 transformation matrices, 
then, 

k' . R = (S R 
* 

= ^ (§ 7 ^ 2 ;) 

= (s’^ 

where S“^ = S’ is necessarily some other member of the set 
of transformation matrices since these form a group* 

Hence it is required to simply look at the way the 

symmetry operations take terms of the kind (r + ) to 

those corresponding to different a s and a s. This depends 

on the way x and x i are defined. Once these are known, the 
oc cx 

appropriate numbers are put into the iEq. (A. 2) corresponding 
to each of the IB Zones and summation carried out. 

A simplification is possible by noticing that changing 
Ic to -k takes H(k) to its complex conjxigate (Ref. Fig. A. 3). 
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Since, 

H0 = E0 - 

implies 

H*0*= E0* 

E being real, the eigenvectors of H*(k) or H(-k) are complex 
conjugates of those- of H(lc) at the same Hi. Therefore, the 
integrands in the BZ integrations for the complex conjugates 
under this transformations k -*• — k. Twice the real part may 
thus be retained and integrations carried out over carefully 
chosen 24 IB Zones., 

From Eq. .(A.3) or from the explicit fora in Eq, (3.32) 
and the s defined in Sec. 3. 2. 2 » the six permutations of 
xyz amount to the following rearrangements of the bond indices: 


xyz 

1234 ‘ 

zxy 

1342 

yxz 

1243 

yzx 

1423 

xzy 

1324 

zyx 

1432 


Also, 


-xyz 

4231 

x-y z 

3214 

X y-z 

2134 


where x, y, z symbolise the components of the k-vectors. Six ' 
permutations of these are again obtained by keeping the 
first index fixed and permuting the other three as above. 
Those corresponding to (-x^ -y, z), (-x, y, -z), (x,-y,-z) and 
(-x, -y, -z) and permutations of each of ■ttiese are clearly 
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olitained by changing the overall signs are taken csn*e of by 
retaining twice the real part of the integrands, as mentioned 
above. 

The required Green's functions for the CCPA calcula- 
tion are. 

Gqq : DLagonal element? r = 0, a = a* = 1, say 
A 

Gqi * Across adjacent bonds meeting at an anion; 
r=0, a=l, a'=2 

Q 

Gqi * Across adjacent bonds meeting at a cation; 

r = r^^, a = a' = 1, 

Gq 2 i Across second nearest neighbour non-parallel bond^ 
r = a = 1, a’= 2. 

where r^^ is the fee sublattice nearest neighbour displacement 
vector, and required parameters in the exponent in the integrands 

r 

are indicated. For Gq^ the exponent is the same as that for 
A 

Gq^ except that a and a are interchanged. This distinguishes 
these. Explicitly stating, if the integrand is looked at as a 
product of two complex factors and C 2 = 0(Sna)0*(Sm*)j 

then those corresponding to Gq^ and Gq^^ have distinct terms like 
Cj^C^ and C^C 2 respectively. 
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Fig. A.l: A tetrahedral part of 
Brillouin Zone. 



Fig. A. 5s Brillouin Zone of an fee 
lattice showing some 
symmetry points. 


Fig. A. 2; Division of a tetralie- 
dron into thin 
tetrahedra (TT). 



Fig. A. 4 8 Irreducible part of 
the Brillouin Zone 
(IBZ) of an fee 
lattxce 0 



Fig. The Irreducible part of the Brillouin Zone of the 

fee lattice and its division into three tetrahedra. 
Coordinates shown are in the units of u/lSa, vdiere 
a is the edge of the direct fee lattice cubic 
unit cell. 
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INVERSE LAPLACE TRAl^SFORMS BY CCMPLEX CONTOUR IIOTEGRATION 
The required inverse Laplace transform is, 

+ ±ca 

a'(y) = / 6^"-® ds 

- ioo 

Considering the complex integral over the closed curve of 
Fig. B.l, we obtain, 

a'(y) =-!-/“ Sin ( 3 “ Sinw) ds 

Since ^ < a < 1, for any finite y, however small, 
exp (-ys) eventually (large s) decreases more rapidly than 
exp (-s°^ Cos Tta) increases. These form the envelope of the 
oscillations due to the factor Sin (s°^ Sin Tta). The distance 
between two successive zeros of the integrand decreases as 
s increases, making the oscillations very rapid for large s. 
These zeros can be evenly spaced by the following substitution. 

Sin Tta = t, when 
. 1 

Q (y) s / exp [-t Cos Tta 

Tt a (Sin 

Sin t dt 

(Sin 

^ * 
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= 2 i„ (y) 

TO (sin 

where, 

I^(y) = exp [-t cos TO 2LL“ ] t(Ve)-l 

”* (Sin f-)^“ 

Sin t dt 

These integrals are calculated numerically by the 
Gaussian integration technique. Fig, 5.2 shows plots of 
\innormalised probability functions Q'(y)’s for various E's 
near the band edge. These fall off to lO""^ of their peak 
values at y ~ 15.0. Hence it is sufficient to carry out the 
integration in a large enough finite interval. For the 
numbers used for the plot in Fig, 5.3 the interval used is 
(0, 50.0), The peaks become sharper away from the band edge, 
and shift gradually to 1.0 as the mobility edge is approached. 



Fig, B.ls Contour of integration 
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ccmditibna (i) a, » i, Oj w 0; (h) Oj » 6, Oj *■ 1. The reeardcm (3) may be vtewed 
aa an infinite set of linear homc^aieona equations in mdcnowns ■{(!«)> Tcet wSS' 
ecHudlstmioy the deterannant 
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The roots of this traneandenial equation are the eigenTaltws. This mode dt writing 
the oharaoteristio equation in the form of an infinite determhuuit m known from 
the eigenvalue problem associated with Hill’s equation (Whittaker & Watson 
iqay). Its apphoations in phyidoai problems, however, had a Bmited sucoms (Kemer 
1951). In the tyidcat case of the anharmonic oscifiatcnr (Biswas d ol. X973) the 
towMt e%ht dgenvaloee were thus obtainable until the numerical errors became 
too sevmre. Tim reason for this is the use of a fixed scale expMxskm (like (t) with 
at M constant) in these works. A fixed scale expansbn is snitalde in a small regime 
of values of (a, A ) where the scale happens to be close to the appropriate vahie given 
hy (2). Outside tlds r^me thn basis becomes unfavourable for expanrion. To fim! 
, the seros o£d{E}ii may be noted that the (hitemdnmtt obtahied 1 ^ mnitting 

all rows mtd offiumns b^cmd the element, he expanded into a fi»B>term 

recundmi between alternate determinants of lower orders: 


With d* *■ 1 normalization, the determinants up to luiy even order (for odd OTders, 
d* «« 0, d| “ i) may be computed successively from the rocurrion ( 5 ). Tl» cfoter- 
lainant d^miE) is a polynomial in If dT degree «i. In the limit of m-i-oo this ^iims 
4(E% The zeros of d«,(lf) are stable for large m (shown in the next section) and may 
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be obtained numericaily by Newton'a method wyeh requires both and 

4m{E). The latter are obtained by differentiating (6) with respect to ff, 

and computing recursively therefrom. The recurstons (5) and (6) used for corapoting 
the eigenvalue are nummcally stable for large m. The accuracy of the eomputa^ns 
is therefore limited only hy the prech^on of the arithmetic used {18 signifimuit 
%ures here). Owing to the quadratic convergence of Newton’s meth<^ it k possible 
to refine a rather crude initial estimate for an tfigurvalue (say, within a small 
percentage) to a 15 figure accuracy in a very short time cm Uie computer. The time 
cm IBM 7044 for the first 100 eigmivalues is a total of two ndnutes. Fch* the 10000th 
eigenvalue the time required is about three minutes. The initial estimates used 
for higher ei^nvalues are the W.K.B. values. For lower eigenvalues a suffic^ntly 
large carder determinant is computed recursively from (5) at various S points. 
Opponte signs of djj for two ne%hbouring E values indicatoi that an mgenvaloe is 
eroBsed which {oovideB a sufficiently acouraie estimate fmr Newton’s method. Hie 


TABUE 1. EtQKBrVALOia O? THB QtrASnO CaCXlJLATOB {B m as*) Ain> 
*H* ABKABMOKIO <}8C1IXAT03 (J? « p*+a:*+Aa:*) FOR A « 1 


quantum 

quartio osoillktor 

snhanooiiie oscillator 

mimfaer, n 

eigenvdueB 

mgenvshm 

0 

1.000 202 000 484 IS 

1.892 351 641 630 29 

1 

2.782 078 028-801 40 

4.648 812 764 212 08 

2 

7.4S5 087 937 980 74 

‘ 8.085 049 957 759 31 

9 

11.044 745 511 378 2 

13.156 803 898 049 8 

4 

10.201 826 018 650 2 

18.057 557 430 303 3 

5 

21.288 372 218 230 0 

23.297 441 451 223 2 

« 

20.528 471 183 082 5 

28.835 338 459 504 2 

t 

32.088 587 710 808 8 

34.040 840 321 111 3 

8 

37.823 001 027 034 0 

40.690 380 082 106 4 

9 

43.f»l 158 087 288 7 

40.865 009 505 676 5 

10 

50.256 264 516 682 8 

53.U8 102 139 685 3 

11 

50.734 214 055 m 0 

60.128 522 858 157 8 

12 

53.402 048 880 718 8 

80.985 030 001 247 2 

12 

70.252 394 028 010 0 

74.035 874 358 102 5 

14 

77.273 200 481 884 0 

81.243 505 060 707 2 

15 

24.457 460 274 942 0 

88.010 348 800 798 2 

18 

81.788 000 808 981 2 

88.128 842 045 234 1 

It 

98.288 800 000 4^ 3 . 

103.795 300 322 273 

18 

100.823 307 381 738 

111.601 815 046 173 

12 

114.080 817 384 985 

118.644 170 733 050 

20 

122.004 638 000 808 

127.017 777 795 355 

2t 

180.042 008 748 030 

185.010 417 826 810 

22 

138.306 147 811 395 

144.142 195 288 390 

22 

147.000 121 26? 004 

162.586 504 805 574 

24 

150.483 502 268 082 

101.144 880 694 513 



Tki anharmmic osciUator 


mi 


TABZiB 1 (coni.) 



quartia oscillator 

aBhansoiBo osotlktlKHr 


aigCRVeluai) 

aigmvalttes 

m 

104.012 043 922 SOS 

100.817 528 001 m 

m 

172.842 711 002 845 

178.600 I«t 800 070 

27 

1S1.S8S 609 too 703 

187.490 ta 692 980 

22 

100.220 238 902 493 

100.488 102 010 221 

St 

mi79 918 883 747 

»>51S32 848 004 

30 

200.232 830 000 144 

214.779 083 840 177 

31 

217.884 201 674 103 

224.074 947 800 

32 

223.833 803 481 138 

288.400 087 470 376 

33 

238.076 280 SO! 000 

242.051 164 081 147 

84 

248.416 303 701 036 

282.628 200 001 4» 

3S 

284.946 107 070 7m 

902.iW 767 400 820 

30 

£04.868 017 814 400 

271.051 8M OSO 007 

87 

274.273 m 030 041 

S81.7M 072 023 830 

3S 

284.068 500 58! 401 

891.723 503 051 018 

30 

203.048 488 260 000 

801.736 843 Ui 187 

40 

303.012 000 348 384 

Sii.831 818 m 701 

41 

SIS.MS 080 183 070 

812.008 000 744 848 

42 

324.008 020 m 133 

832.284 003 530 091 

43 

334.201 782 m 402 

S4SUS90 878 832 547 

44 

844.877 020 801 886 

388.012 000 233 780 

40 

884.080 688 600 308 

363401 804 008 470 

40 

305.378 444 407 008 

874.000 879 «}0 008 

47 

878.802 863 004 033 

384.705 074 675 721 

43 

380.480 330 OSO 817 

308.415 048 465 203 

40 

307.141 320 780 074 

400.200 007 443 tU 

00 

407.874 803 284 438 

417.080 203 284 848 

too 

1020.980 092 108 37 

1086444 183 138 91 

1000 

2! 866.202 118 137 7 

31032.783 710 000 0 

10000 

470700.204 427 028 

471 103.777 700 800 


•igtt diange of a sufficientty large (»der lieteminMit between two aei^dKmiiag E 
valuee may be used saoceedvety to yield increaeingly aooorate o^ier and lowm 
bounds for the eigenvalues. As a method for the oomputa^km of dgeovdues tha 
is much slower than Newton’s metiiod. However, mi eigenvalue computed hy 
Newton’s method may be checked ahd in the procen be upper and lower bounded 
l^.the sign check desedbed above. The eigenvalues computed in tUs work were 
evaluated to 16 significant figures, upper and lower bounded in the 16th fig^, mid 
Uienroundedoff to IS figures for the tables 1 and 2. Table 1 lists smne of ottf results 
for the anharmmiio oscillator (N «» |>*+aj*+Aa:*) said the j^uarttc oscillator ; 

Each ' % 'nvaJue in the tables has been recheolwd by udng sevmul whies of a in 
theran,^«. «-|+(t.2 to 1.4)iiW, (7) 
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in which the constant within brackets {t .2 to 1.4) is set empirically. The effect of 
using different os from (7) in the computation of eigenvalues is that the stabilization 
of the initial estimate for an eigenvalue to 16 figures occurs at a slightly different 
point in the recurmon. It emphasizes the point that though the scaling must be in 
the correct range for an effective expanriim there is no exact scaling neceMKry. 
Table 2 demonstrate the uniform applicability of this method for all value of * 
and A. VVe focus on the r^sime of extreme value oi (n, A). Hioe ei at. (1975, 1976) 
give different formulations of the eigenvalue proWtem in each of thee regime, 
which cannot be extmtded into the other r^me beeue of the ‘boandwy layer* 
in between. The eigenvalue omreponding to n » 1000,A« 10~* is in the ‘boundary 
layer * and also in the awkward h%h n, low A r^j^me (Hioe et al. 1976). In the present 
approach this eigenvalue is obtained by the eme technique and with the same 
accuracy as any other, highlighting the scope of computation with an a^idPOiaiately 
scaled bads. 


Tjcbxx 2. EiosinrAi.tn» cnr thb anbabmomio <»G nAA.TOs {H » 
m zuiomis 07 bxtrsus values ov (a, A) ■ 

A 0.0001 40000 

n 

0 1.000 074 686 880 20 SS.274 408 188 780 8 

(near hsrmonki r^iiine) (near quartk t^me) 

1000 2134.242 648 232 21 747783.421 602 834 

(boundary layw) (iMarqttartic/W.K3. 

rtg^noe) 

8. Stabilitv ov zb&ob ov 

The reoursive evaluation of the determinants dJiE) and the staUUty of thdr 
zeros for large m will now be considered in some detail. It may. be noted that die 
recursion (6) is obtained from the recorsicm (3) by replacing with and ohuiging 
the sign of every alternate term. This prescription is vdid when the ooefildent of 
the highest order term in the {o,*} recursion is set unity (by j^operly divuiix^ if 

^,.(-1^ m.OXt...:] , 

The asymptotic behaviour of the solution (rf SchrSdinger equation for the an- 
harmordo oscillator approaches exp ( ± which gives 

for the coeffioiehts in the series solution { 1 ). It follows immediately from the reladcms 

OT-*-CO. (10) 

The a^mptotio relation (fO) implies a sequence of decreasing determinmits beyond 
a sufficiently large in. This zerdng for large in must be isolated from the ^t«- 
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laiaatiwi of the eigonvahiM which are the aeroBof for values of U. It ispossible 

to achieve this quite simpiy by multiplying the reoursimi (8) by a laige number 
wbenevOT, while applying I^ewtcm’s method, the reoundvely computed determinanta 
become too sm^t in magnitude. This zenorms^ization amounla to starthig the 
recurmon with o h^;hmr value of the arbitrary constant Mare generally, other 
recursively connected sequences of determinant may be define such that 
the Eeros of *ttd ate oomnum but JlJiB) may be givm any deMred 
asymptotic behaviour for liwge m. This is done Iqr multiplying the recundon (3) 
by a function of m, say /(m). The ocmreBpcmding infinite <fot«rminant A(E) and its 
various ordcae truncations are related to the respective quwitilieB for A{E} Iqr the 
relation 

4(fir)-/(iti-2)/(«-41.../(0 or (il) 

where tho right hand aide contains /(O) (or/(i}) for the even ((W odd) elgraividoes. 
Oleorly, can be given any asymptotic behaviour for !azp m by pn^ly 
choosing the fonction/(in). Sina>/{m) is independent of £ 1^ d^bitioot tim exm 
of ^miE) and AJ(B\ are common and they are equally suited for the emnputation 
of the eigenvalues. The ‘renormalisation suggested above is a ^e<dal case of 
this multiplication’^ in which all rows of A{B) are left intact except one udddi is 
multiplied by a large number. 

To see the stability of the zeros ot A^(B) as m-^oo we first con^iar the same 
problem for the harmonic osdllator which is exactly soluhle. In the case (dt 
harmonic oscillator the zeros of the characteristic polynomial At^iE) are real mmI 
the ratio of successive polynomials is 

(am-t-i-J) 

4(1) ” (ffl+iy(ia+ir 

Tho zeros of A^^.f,(B) consist of all tbs serot of ^ ■» 2«+ 1. 

Tiiu® tho eigenvalueB resulting from the solution A^(E ) » 0 are reproduced 
exesdlyl^ the solutions of " 0 for all «i. This is characteristic of miexaetiy 

soluble problem. For anharmonio oscillator the zeros of successive mder pdy- 
Koaaisfe SEO different. For a given F and suffidently laige « however, 

; 

d«+,{F)/d*-,(l?) ~ -A/w*. m-^co. 

This implies that the langeat wder t«m in the ratio A^,(E)/AMiS) ia ittd^elMhMtt 
of E for The successive polynomials as functions of E therrfore ^er by a 

mcltipUcotive constant (depending on w») for large ui. Hence, the zeros 
Btabilke for large m. 
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4. ElOBHFrrSOTIOKS AKT> MOHSHTS 
The recursion ( 3 ) falls off for large m according as «m+i/a»,_4 ~ 

This ^sures a rapidly decreasing (in magnitude) set of coefficients {o„} after a 
sufficiently large value of *«. When M is set equal to a computed eigenvalue in 
recursion ( 3 ) the resulting coefficients provide a very convmient representa- 
tion for the corresponding eigenfunction through the expansion (1). To test now 
well the computed eigenihnctions satisfy the Sohrodinger equation we compare the 
two sides of the test equality, Hi^(x)l^{x ) » E, at various points *. For the first 
ten eigenfunctions which we have computed the test equality is satisfied to 13-14 
mgnificant figures from x « 0 to points well outside the cktsacal r^on. For jij co 
the right hand side of (1) must go to +oo or —oo unless the value of B used in the 
computation of the rolution is an exact eigenvalue (Tltohmarsh 1961). The blow-up 
occurs in the non-clasrica! region and shifts to largex asE api»roaches an eigenvalue. 
However, this large x behaviour does not affect the computation of accurate eig«n- 
funotions significantly since well before the blow-up the computed eigenfunction 
reaches extremely small values in the non-classica} region. The part of the computed 
eigenfunction where the blow-up occurs may therefore be replaced by xero without 
losing much information. In the typical case of the lOth eigenfunction of the 
pure quartic oscillator the test equality » Eiq is satisfied to at least 13 

significant figures in the entire classical region (x » ± At a pmnt x » 1 .SBff the 
test equality is still satisfied to 10 significant figures, where the value (ff the computed 
eigenfunction ^^^x « l.SFfo) is O(I0-w) relative to « 0) » 1. The tranrition 
moments between the enharmonic oscillator energy eigenstates, g^ven by the matrix 
elements may now be obtained from the computed eigenfunctions. It is 

shown in a recent brief report (Banerjoe 1977) that thew matrix elements satisfy 
an exact liireor recurrence relation with respect to the index k : 

4*(fc - 1) r ji> 

-bi!(h-l){h-2)(t-3)<il**'-*|7>, . (12) 

where the potenticd V{x) could in general be any polynomial function. P<w the 
anharmonio osdilator therefore the lowest moment of any particular traiuntion is 
sufficient for the computation of all higher moments (of that transition) recur- 
sively, nwing the eigenvalues obtained above. The computation of the matrix 
elements <»|*|j> and <*(**|i> is done as follows (i) each eigenstate is defined by 
including as many terms in the expansion (1) as are needed to obiidn the cone* 
spending eigenvalue to 16 significant figures; (ii) as described earlier, tl» computed 
eigenfunction falls off to very small magnitudes in the non-classical region befine 
finally increasing for x-»-±co. The range of integration for the ©valuation of 
the matrix elements is therefore truncated at a point at which the computed 
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Tarlowa loealisatioa crttf«riw %!si*t ls*r« 4* tfes 

p«ni.h*vu rwc'jt’pnf'aijr ebown tlse «5L»iwo« at *l®«troR 

i.a di®<»3r4®rtHi but fall to «faow Xh% 4*?ort»sit f*atu 7 ’ 4 !., 

■vlu* tfle trisfiaitloa to Blocb t^pe aoa-iWBSiaiiBafei* atatfi* **sroii 9 lAs 
woMUtf «<g#. la th» tre&tamt of Abou iSsaeza aobtli^ 

«So«s not approach baa^ »4g® 1%, ‘ta t^« <ll«rttr 0, 

Thr crttfrla auggestad by AM^m^nk2) mA ishtiu-, rigotmw 

4n principle are, ssathmatie ally intractable kwplug m® rigour. 

Brnd tic simplifying aasuroptioas ar« therefora W® bar® trl«4 to 

look at toe radius Of locallsatioa la tb© localiaM asrwgy r*f, 5 loaa. 

W® t,*»» aa th« delooalis&tion criterion Itaelf, th« prctjabUiatic 
dl’f'‘8Tf:encij of the loealisation dcusain* ?or tb* wsak disorder 
■|^» mobility edge Zq Ib thus located within the bamS (this htus been 
raptart««J by X’jai&i* et.J:.v4J) «nd, 1® obtedaed aa a fujjctioa of the die- 
order paraKieter. Sho analysis in the me® lAQa disorder is larg®** 
leedlf^ to Anderson traaeitioa is possibl®. yurther» the behsrlour of 
the «i*e of the locallsatioa dooaln and the way it diTerge* 

e« fc }$c i.s looked «t» The localisation length is a panaaeter that 
rtBttre the expressions for yayioias physical properties eg. phonoa 
«na4at«hl hopping eonduotlonw) and her*® the aotiwatioa for the mXm- 
lation for the same. 


II , rc:?:.naATicii 


Con:}i4er the Anderson !l!i»iltOKd.«n 
1 " 


it) 


aM repx’wseatation Of tb» r«8olr«at (ZI-hJ . If «n electron 

le p'U'r't^S at site % at ttae t » 0, the probability that it is there at 

?„ • W ••• isj 

if C,,,,.;*) * [a- o “ Sc(»)r'‘ m^.'oxsrnuoa 

beixjg the self energy fUtcUiia, then 

1 dl ('~ fj) * 




A (l« .'-sCi) f> O'j, for *rt»»ad«iS Aif-nv/M *.i,4 

V « , whac for :oc«Iia’xi atiite* a*cl «.rc ^ 
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^o/ 1 , 

1.. .*# tons ,i4e&t:Ty 



-4 1 .-f 


1 

i0 


}** 0 ^i/hi tor « l 0 <?«l 4 «.d «T*p t 

I*«t F I4} ib« •■S iHift'ributloai m^v S^t 1, 

itfhmm 2j ^ d%/ dK/ thm F{3t» 3j, Ita Four! »>r-Wpiae« %3tm»iotm^. 

fouM H fwliowtai; non-U -Motr iatafiwii 

'*« " *'*' I 

f* f'lrat attnijpt a aolutlott of th« v^rp# 

F(fe,S)a 

nhrm d«aot«« « aot of parMotor* dowisriljlaf «»# pwPisJiiUtty tlm- 
«ity of »Mei)i 1/i iMepoadont of tbo psofeateai^ doiwlily of 1 } *04 
TOt id 00 rmt*.* &• 'r8*diit»l«« Ij *sad % isjte* horwft &9pm4.mi» a# 
abw® fowl thu» doao mt iMpXj rtvitloml tis&iipmimm of Bi iwi 1«, 
l»t i« 4m appjpoaslWitlou noao-tho-l#**. la partiowliMPt 


( 4 ) 


( 5 ) 


f (h,&3 •» 

Siletaf p(4) to iNi a CiimdhEf diotxilmtiostf w$ htcrm 

|(»l) 

^•404 P(i) Mi'tiiiiyii%' 1Ai« foULowiail oquation 

% «x»a % bolng glToa ' 


•do#' 


(«) 


(f) 


^4*. + ^ '^t' * 






t * 


..V (f) 


F(jl) ommot too nol-rod for aaaJ^dlesJlly tem (s) boomuio of th« 
tap (-SF^A^) ftetor- But la th* otak dioordor limt It (f+ai)«,i 
It adaltt tht foUottag tc^utlmi 


(S> * e*** 4a4» 


oB-tf* 


k¥‘ 




^ a 

4 


tn (k) 




!tl' 


Ci) 2? 


... (?o) 

... (ti) 

» •• (if) 


IftCr)# ttJt prokmixtf 4i«ti4but3.0B of fist lootUattlcw is tis« 

ofetsUatd to tiwi StplAOt iwtrt* of FJ^Cs). ?&««« »■*•• eyrrtt of tii» 
»i»«o ia Htf. 1* »i»«t t»# *»«M# t,'loo« 4«il, Ch,} 2.* f««« to 
ft tiMUf dtftat tAt armetcfc iowjUftotioa *» 


H13 


,# • # 


m} 


^ V 1, «ai f(s) tppscwdM* • i <fuaBfl«iBa 

<u'>-’-*|( “■ *'] { 

®i® ot a«itt' Mgr b* •*«b 1» b* liaatar «• vtnilA It* 

«q^et«(i frem «. R4»*i;i fliilA 'ystoxgr. Sh* €S& 1« alda to 8U«b aa 
«i®itsoajmattoa» ®i« BsiBaifyl*^ r»»«lt« ar« niioiai in Fie. 2. 
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But as «• increases this compromise scaling becomes unfavourable for taoxt wid 
more eigenvalues. The way out is to compute each eigenvalue separatrfy by using 
en appropriately scaled set of basis functions. This is intractably labtmous in a 
variational scheme. In our method the use of an appropriately scaled basis mwely 
requires that a proper value of a obtrined from the formula ( 7 ) be used in the 
recurcions (3), (5) and ( 6 ). Since each sigenvslue is computed individually tlmie 
are no carryover errors. 

(v) Computation with a larger basis Is very simpfy done in thhi inetl«)d by 
continuing the recursions {3), (5) aand (3) for increasing m. In contrast, a variatiomd 
calculation with a larger basis requires integration and the subsequent diagcmaliza* 
tion of a large matrix which beyond a certain size is intractable. For instance, the 
lOGWth ©igenvaluo of tho anharmonic oscillator stabilizes to a ifi Hgurs accuracy 
(in three minutes <m IBM 7044) at a point in the recursions which ccareiqmnds to 
the use of nearly 1 7 500 terms in the expansion { 1 ). A variational calculation of this 
size is inconceivable. 

(vi) The eigenvalue problem of the general anharmonio osciBatar 

j5f la jpS + aj® H- 

^ w 2,3,4,... is solvable precinely along the came lines. The generidized scalii^ 
formula is 

a(n. A) « H (S') 

Thus the following regimes may be distinguished; 

j is the near harmonic rdgime, 

Id the near pure anharmonio regime, 

jjt»-i); 0 (+«^i«/i+j> « I }g the boundaiy laycw between the above two 

regimes. 

Tho importance of tho combination »^~‘A in detwmining the above regimes wm 
recognized by Hioe et al. { 1976 ) on essentially empirical grounds. It is in the present 
work that the combination rt^~*A is shown to determine the charactermtic scaling 
in a given regime of (n,A) tlrrough the relation ( 2 '). This observation leads to the 
construction of a scale adapted basis and to a uniform treatment of the problem in 
all regimes of («,A). The detailed results on the general enharmonic oscillator 
problem mil be reported later. 

Computational assistance from Mr M. Ramesh is gratefully acknowledged. 
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( 7 ) 

occur in the 7x7 matrix G • Interactions between these 

bonds specific to each of the constituents (GaAs and InAs in 

case of alloys) are used to determine the elements 

Of the 7x7 matrix H for each configxiration and then the 

M 

7x7 matrix G is determined from Eq. (3.28). Finally, a 
weighted average of the (1,1) element of these (since the 

central bond is numbered 1 ) yields the configurationally 
averaged diagonal Green's function sought, 

(e) Finally, a cluster of six bonds, forming the characteristic 

ring in diamond structures, is considered. The configurations 

and their respective statistical probabilities of occurrence are 

shown in Fig. 3.7(e). Again, the 6x6 matrices for 

and G^^^ are written down and G^ determined from Eq. (3,28). 

M 

An average over the trace of the weighted averaged of the G 
corresponding to each configuration of the cluster is taken. 

This is the diagonal Green's function sought, averaged over 
configurations of the ring. 

The effect of such ring clusters is important. In amor- 
phous alloys, when the orbitals are distorted, we often get five 
and other odd membered rings. The effect of such rings on the 
density of states is a matter of much interest. It is important, 
therefore, to first study the effect of rings in the crystalline 
materials and alloys. 




Fig. Mil (i) -1“ md (11) -TM o“’* Wr 
0^),5 As uslag the fitted 

paraaiettrs* 
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This effect of non—oi’thogonality has been argued, away b 3 ^ 

Chen and Sher but their arguments are weak and difficult 
to QUStif^r, 

In this chapter, an ab initio chemical pseudopotential 
scheme originally introduced by Anderson (lb 6 9) and later 
developed by Bullett (1975) is used with the aim of over- 
coming the above shortcomings. Herman-Skillman tables of self 
consistent Hartree-Fock atomic orbitals and potentials are 
employed as the starting point to determine the pseudo- 
Hamiltonian matrix elements characterising the III-V semicon- 
ductors under consideration. While no claim of extra-ordinary 
accuracy of the parameters so obtained is made, the procedure 
is superior to parametric fitting since it does not depend 
on experimental results. Further, non-orthogonality of the 
basis set is present naturally in the scheme and no approxi- 
mation is made on its account. This is where the main thrust 
of the scheme lies. 

The procedure is described in the following (Sec.A.2). 
It was pioneered by Austin (1962) and Anderson (1969) and 
extended for application to the cases of homopoiar carbon and 
silicon in the bond basis by Bullett (1975) » The chemical 
pseudopotential scheme has been further generalised in this 
work for application to heteropolar covalently bonded zinc- 
blende alloys. In Sec. 4,3 this is applied to the cases of 
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calculation on benzene in the site orbital basis. This consists 
of defining the Wigner-Seitz cells corresponding to each bond, 
considering these to be situated at the centres of the lines 
connecting the neighbouring bonding sites. Operatiorially, this 
is done by checking that each point in the Wigner-Seitz cell 
corresponding to a bond is not closer to any other bond centre 
than its own. The 'Wigner trick' involves assuming that an 
electron at any point sees the potential due to the cell in 
which the point is situated only, v;hile the potentials due to 
all others is screened by the Fermi sea. This is equivalent 
to dividing the crystal into Wigner-Seitz cells, and the 
potential in each of these is only its own bond potential. 

Under this approximation, the pseudo-Hamiltonian matrix elements 
are like 


<a' o' |ao> = / ?>, 

o’ ' 


a c 


<Vo' 




ab 


d^ 


= StAsT J V^a> ^a ■ • 

"a'o' 

+ iVc'-"o’) «*ab3 

the integration being over the Wigner-Seitz cell Wg^,^, corres- 
ponding to the bond orbital la'c'> only. 
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four alloyed cation sites respectively, the seven bond clus- 
ters Gy and the six bond ring clusters Cg are embedded in 
the 1-CPA medium as an approximation to the corresponding 
n-CPA self-consistent medim as before* Density of states 
per bond (in units of .corresponding to each configuration 

of each cluster and the configurationally averaged ones are 
calculated and graphical results for GaQ^ 2 . ^^0.9 pre- 

sented in Fig, 4,5 and in Fig. 4,8 to Ga^ ^ ^^0,5 ^ ^ times 
the density of states per bond corresponding to the VGA and the 
1-CPA are shown on a larger scale in Figs. 4, 6 and 4.9 for the 


same alloys. In Figs. 4,7 and 4.10 again, (a) -]jn 

(b) -Im«G»p , (c) -Im«G»p , (d) -In«G»p and 
C4 Gy 

(e) «G»„ , averaged over configurations of clusters C-,, C^, 
C^i , Cy, Cg (See Fig. 3.7(a) to (e)) respectively are shown 
in a larger scale for clarity. 


The comments regarding these resiilts are exactly analo- 
gous as those for the fitted parameters. The effect of clus- 
ters is the same as before. In the x = 0.1 case there is an 
enhancement of the upper lying states, arising mainly out of 
isolated clusters immersed in the medium. There is now also 
some effect in the s-like part - again an enhancement of the 
upper like peak. The 50 - 50 case shows less effect, but with 
perceptible filling in of the upper lying states near the 
upper band edge. 
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of disorder «», This insnediately suggests the following 
identification of the size of the localised state, 

Lo<E) = (1 - ) (5.3) 

The identification is justified if we also note that the 

^^0 ■*! 

residue of Gqq at its pole is (1 - -gg) • 


For states localised over L sites, using the proba- 
bility interpretation, this is also ~ 1/L. Hence a justifi- 
cation Unfortunately, as we shall see, the proba- 

bility density of Lq(E) has a long tail leading to diver- 
gence of <Lq(E)>. However, <10”^ (£)> exists and it is 
more profitable to directly look at this quantity and 
definetie averaged size of the localised state as 
Lq(E) = <Lo“^ (E)>"^, 

5,3 The Probability Density : 

We shall study the simple model of diagonal disorder 
with non-random, nearest-neighbour overlap. Within the 
Anderson tight-binding Hamiltonian, 


H * £ e. P. + E V T. . 

i ^ ^ i,j 


V. . = V 6 (r . - r . - x) 

1 J U 


where x is a lattice vector, and form a set of statisti- 

cally independent, identically distributed variables with 
probability density p(e). 
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Qo(y) S 6[y - Sec 

when, 

Sec + 1} f5.20) 

We note that near E = E^, <Lq’’”>“^ diverges lirearly. 

This is the familiar result we would expect from a mean 
field theory , The CTA is akin to such an approximate 
theory. We cannot conpare this with Anderson* s prediction 
of a 3/5 power-law divergence. His argument was based on the 
study of on a real lattice, and the 3/5 power law 

near E depended on the statistics of closed, ncn-reptcting 
paths in the lattice. Such a result cannot be expected ’..ithin 
our CTA. Departure from the linear law is evident away frem E^,. 



Fig. 1. Averaged value of the inverse size of localisation do- 
main in die iocalised.reginie for a veiy weakly disordered sys- 
tem. Site energies are Lorenzians with very narrow width 
t/u«0.01. 
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of the square face. The IBZ is then l/48th of the BZ since 
there are eight hexagonal faces and six square faces. The 
tetrahedrons FLAK, riUA, ruXA into vdiich the IBZ is divided 
are of volumes 32, 32, 32 2/3 respectively in units of (u/lSa)^, 
For the calculation, each tetrahedron is subdivided into 64 TT's . 

I 

since finer division does not reveal any additional features 
in the Green's functions evaluated. Energy eigenvalues have 
been determined at five points besides r, the tip of each TT, 
spaced evenly along axes of each TT and 6-point Lagrange's j 

I 

interpolation formula ! 

j 

f(x + ph) = - Efe?rl^ig-2) (p-?) 

i 

. f + j-it) .(^2l ! 

12 ^o ^ 12 ! 

i 

is used where f^^, i = -2 to 3 are the known values of the func- i 
tion f(x + ih) for i = -2 to 3» i.e. at intervals of h around x. 
The values of such that E = have been obtained by Newton 
Raphson's method. 


For the imaginary part of the diagonal element of the 
Green's function in the coordinate space bond bases set jr a>, 
f(k) s= 1 and the above procedure is complete. In general, 

ik. 


I 5 

k.n 


"tta’ (A. 2) 


E- 


<ra |G jr' aV 
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The Enharmonic (^dilator 
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iiAowiSite <^gwiv&Itteti and eigenfanct!<»u cf tha anhwmonie ORdlktor 
{H h>ii) and the qear^o osdllstor (H aw 

obtained in ali rhinos o£ the quantum number n and ttwiw monifltty 
cmastant A. l!ranaition nsomente of compsrabfe aooura«y jure obtained 
for the ({Qartio osdilatw. 

The nwthod, applicable quite generally fcr ea^Yslae proldems, k 
non-perturbathre and InvolTea the ueo cf an apjiropriately aoaled batis 
for the determination of each mgenvalue. The ap|»oi^te soatk^ formula 
for a ^en rtpim of (n, A) is ocnetruoted from the osdllation properties 
of the eig^nmnotions. More gmieral anharmonio c^oillaton are also 
discussed. 


1 . iHTBOnCOtlOS 

The quantum mechanical enharmonic oadllator dmcribed the Hamihonkm 
jyr wjjs 4.3^4. Ao^, where p«>id/da; and A > O ktlmanharmonidty constant, Ime 
been the subject of very geoeral interest. The inves^kons fimm tlm pednts oi 
' view of moieoular physios and field theory are refemd to in the works of Reid ( 1970I 
and Bender & Wu (1976) re^nrathrely. A general account may be found in Hioe 
A Montroll (1975) and Hioe, MaoMillen A Montroll (1976). 

The perturbation expansion for the dgenvalues ii 7 «(A) in powers of the 
harmonicity constant A is not convergent but a^mptotio (Bendw A Wu 1969; 
Simon 1970). Such investii^tiona have been confined to the computatimi oi 
lowest few eigenvalues for various A. Variaticmal cidouIati<ms uidng a hanmmio 
oscillator basis have been more successful Beid (1970) tht» obtahied &st 28 
eigaovalues of the quartlo oscillator «p*+af*) to 12 sigidficant ^ures. For 
higher eigenvalues the variational calcutadcms beeome forbiddingly laborioos. 
Moreover, the variatiohally obtained eigenfunctions are far less accurate tlum the 
corresponding eigenvaiuefl and are therefore nnsuitable for the computatioa of 
matrix elements of operators. Biswas cl ol. (1973) uasd the Hill detenainent method 
to obtain eight eigenvalues dt the anhamionio osdllator for various A in 'the nmge 
0 < A < IOO4 For highcar eigepvaluea or fw higher A the nuR»rii»d em»s in th^ 
work become too severe. The most comprehensive wwk rat this pfdblem it due to 
Hioe ft Montroll (1975) and Hioe et al. (1976); they distinguished two limiting 
regimes of valuee ^ the quantum number n and the smhwrmraticity eimstant A. In 

I m I 



